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A  new  experimental  approach  to  the  determination  of  three  flame 
diagnostic  parameters  is  developed.  The  experimental  application  of 
saturated  atomic  fluorescence  to  the  measurement  of  quantum  efficiency, 
total  number  density,  and  noise  power  is  discussed.  Data  for  quantum 
efficiency  and  total  number  density  are  compared  to  the  scant  litera- 
ture values.  Four  elements  and  three  different  flame  compositions 
are  investigated.  Preliminary  noise  power  spectra  are  discussed  as 
background  for  the  use  of  saturated  atomic  fluorescence  to  measure 
noise  in  the  atomization  and  nebulization  processes. 
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CHAPTER  1 
INTRODUCTION 

Characterization  of  physical  and  chemical  parameters  of  flames 
and  plasmas  has  been  a  goal  of  chemists  (1,  2,  3),  physicists  (4,  5,  6), 
and  combustion  engineers  (7,  8,  9)  for  many  years.  The  goal  of 
physicists,  physical  chemists,  and  combustion  engineers  has  been  to 
understand  on  a  quite  fundamental  basis  the  processes  which  go  on  in 
flames  and  plasmas,  i.e.,  oxidations,  reductions,  free-radical  reac- 
tions, ionizations,  etc.  In  recent  years,  an  even  larger  impetus  has 
been  seen  for  this  research  as  sources  of  economical  energy  have 
dwindled  and  better  combustion  understanding  and  design  have  become 
paramount  in  importance.  Analytical  chemists,  on  the  other  hand,  have 
become  interested  in  flames  and  plasmas  because  of  the  extremely 
important  role  of  these  devices  in  quantitation  of  elemental  species 
in  various  samples.  The  coincidental  research  interests  have  brought 
all  together  to  try  and  provide  a  better  understanding  of  basic 
processes  common  to  both. 

Quantum  Efficiency 

A  quantity  of  great  interest  to  the  analytical  chemist  and  other 
scientists  is  the  quantum  efficiency,  the  measure  of  the  relative 
yield  of  fluorescence  processes  to  absorption  processes.  The  quantum 
efficiency  has  a  strong  influence  on  the  performance  of  normally  (low 
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intensity)  excited  atomic  fluorescence  spectrometry,  which  has  become 
an  important  flame  diagnostic  (3,  10,  11)  and  analytical  tool  (12,  13, 
14)  in  recent  years.  In  addition,  quantum  efficiency  is  an  indication 
of  (rates  of)  deexcitation  processes  in  flames  and  plasmas  which  can 
be  invaluable  for  unraveling  many  basic  spectroscopic  and  fundamental 
properties. 

The  traditional  method  of  measurement  for  this  quantity  (Y)  has 
been  to  measure  atomic  absorption  from  a  continuum  source,  then  to 
measure  atomic  fluorescence  with  the  same  source  positioned  at  90 
degrees  to  the  absorption  axis.  Great  pains  had  to  be  taken  to  ensure 
that  solid  angle  effects,  flame  edge  effects,  etc.  were  eliminated 
(15,  16).  At  best,  this  yielded  a  spatially  unresolved  picture  of 
deexcitation  processes  (i.e.,  absorption  is  a  line  of   sight  tech- 
nique). In  addition,  no  hope  could  be  held  cut  for  temporal  resolu- 
tion. Also,  the  measurement  involved  two  highly  different  electronic 
gains  (those  typical  of  absorption  and  fluorescence  measurements), 
which  must  be  calibrated  against  each  other. 

The  measurement  of  quantum  efficiency  via  saturated  atomic  fluor- 
escence, however,  allows  both  spatial  and  temporal  resolution  under 
the  proper  conditions.  It  is  inherently  easier  to  measure,  as  it 
involves  only  ore  optical  train  for  which  the  solid  angle  must:  be 
found.  However,  the  photomultipl ier  tube  and  optical  system  must 
be  calibrated  via  a  standard  source.  Along  with  the  system  used  for 
this  investigation,  (a  nitrogen-pumped  dye  laser)  comes  the  need  for 
quite  sophisticated  electronics.  Signal  processing  and  cross-checks 
for  problems  become  much  mors  involved. 
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Total  Number  Density 

Total  number  density  has  been  another  sought  after  quantity  in 
flame  and  plasma  spectroscopy.  An  important  figure  of  merit  in  atomic 
spectroscopy  has  been  the  efficiency  by  which  a  particular  atomizer  is 
able  to  convert  solution  or  solids  introduced  into  it  into  the  form  of 
neutral  atoms.  In  most  cases,  this  efficiency  has  been  separated  into 
two  parts: 

1.  efficiency  of  nebul ization-yield  of  fog  droplets  from  a 
volume  of  solution  introduced, 

2.  efficiency  of  atomization-yield  of  neutral  atoms  from  sub- 
microscopic  species,  i.e.,  atoms,  molecules,  ions,  etc. 

The  former  quantity  is  a  very  tedious  one  to  measure,  involving  long 
periods  of  aspiration  of  a  relatively  high  concentration  solution  of 
a  suitable  element,  collection  of  waste  solution  and  washdown  from 
the  nebul ization  chamber,  and  dilution  to  volume.  Then,  the  concen- 
tration of  the  analyte  in  the  diluent  must  be  determined.  Typically 
(for  most  nebulizers  in  atomic  spectroscopy),  this  is  a  low  efficiency 
process.  Therefore,  the  accuracy  is  subject  to  determinate  errors  in 
procedures. 

The  problem  of  determination  of  efficiency  of  atomization  follows 
directly  on  these  results,  unless  a  separate  procedure  is  developed. 
This  is  virtually  impossible  unless  a  vapor  generation  technique  is 
available.  The  efficiency  measurement  must  be  made  via  detenrri nation 
of  the  total  number  density  in  the  flame.  The  classical  method  for 
this  measurement  has  been  to  use  absorption  from  a  continuum  source 
(17,  18),  which  provides  the  desired  quantity  only  after  a  tedious 


calibration  procedure,  including  correction  for  reflectivity  of 
multiple  surfaces  in  the  optical  collection  train.  An  expensive, 
high-resolution  monochromator  is  also  required,  as  the  atomic  absorp- 
tion line  is  quite  narrow.  In  addition,  this  method  does  not  allow 
spatial  or  temporal  resolution. 

Saturation  of  the  pertinent  atomic  transition  provides  a  fluor- 
escence signal  (therefore  a  spatially  resolvable  signal  in  the  usual 
configuration),  that  is  not  a  function  of  flame  conditions  (except  as 
a  function  of  atom  production  capability),  and  is  a  function  only  of 
basic  atomic  properties. 

Noise  Power  Density 

Noise  power  density  measurements  have  been  used  typically  by 
information  transfer  and  electrical  engineers  as  a  means  of  discovering 
the  sources  of  systematic  and  random  noises  in  hopes  of  eliminating 
them  and  being  able  to  transmit  data  with  fewer  errors  at  ever  in- 
creasing rates.  The  techniques  of  noise  measurement  and  theoretical 
treatment  have  been  growing  in  importance  for  analytical  chemists 
(19-28),  as  the  science  grows  in  depth  and  sophistication.  The 
attempts  at  measurement  to  date,  however,  have  been  troubled  by  lack 
of  definite  calibration  and  lack  of  willingness  to  track  down  and 
eliminate  sources  of  instrumental  noise  in  analytical  techniques. 
The  difficulty  of  this  task,  while  formidable,  is  not  insurmountable 
(29).  If  calibration  can  become  complete,  and  the  sophistication  of 
analytical  chemists  can  be  improved  in  areas  which  are  now  subsidiary 
(flow  engineering,  electrical  engineering,  etc.),substantial  improve- 
ments should  be  forthcoming. 


Saturated  atomic  fluorescence  stands  to  provide  a  substantial 
contribution  to  noise  analysis  in  atomic  spectroscopy.  Saturation 
provides  several  features: 

1.  elimination  of  dependence  on  quantum  efficiency  (flame 
conditions) , 

2.  spatial  resolution, 

3.  possible  loss  of  dependence  upon  the  atomization  step 
(if  correct  analyte  is  chosen), 

4.  possible  relative  freedom  from  emission  background  noise 
(if  conditions  are  chosen  carefully),  and 

5.  freedom  from  source  variation. 

Therefore,  saturation  can  further  the  understanding  of  noise 
sources  in  atomic  analytical  spectroscopy.  It  can  provide  separation 
of  noise  of  nebulization  from  noise  of  atomization  (if  such  a  noise 
exists),  allowing  better  understanding,  and  design,  of  nebulization 
systems. 

Reaction  Rates 

Since  analytical  flames  and  plasmas  are  used  at  atmospheric 
pressure  and  are  quite  hot  (2000  K  to  6000  K),  reactions  are  myriad 
and  fast.  In  most  flames  and  plasmas,  some  sort  of  local  thermodynamic 
equilibrium  (LTE)  is  assumed.  In  several  cases,  however,  the  mechanism 
of  excitation  of  atomic  species  is  unclear  (e.g.,  the  induction  coupled 
plasma,  which  is  not  in  LTE).  Engineers  and  physicists  share  with 
chemists  the  need  to  understand  these  and  other  phenomena  kinetically. 
Saturated  fluorescence  (atomic  or  molecular)  can  be  a  very   useful 
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tool  in  this  pursuit,  as  a  means  of  simplifying  kinetic  schemes  by 
swamping.  Reactions  such  as  photoionization,  quenching,  and  inter- 
system  crossing  can  all  be  studied  by  means  of  saturated  fluorescence. 

Tunable  Dye  Lasers 

The  advent  of  the  tunable  dye  laser  (30,  31)  and  the  power  it 
gives  to  observe  atomic  and  molecular  populations  have  given  flame  and 
plasma  spectroscopists  an  invaluable  tool  for  diagnostic  procedures. 
The  general  properties  of  these  and  other  lasers  include: 

1.  directionality  (low  divergence), 

2.  monochromaticity, 

3.  coherence  (spatial  and  temporal),  and 

4.  high  irradiance. 

All  of  these  laser  properties  have  been  useful  for  analytical 
atomic  fluorescence  soectroscopy  since  Fraser  and  Winefordner  (32)  first 
used  a  dye  laser  to  excite  atomic  fluorescence.  Their  work  covered 
nine  elements  in  hydrogen-air  and  acetylene-air  flames.  Limits  of 
detection  obtained  were  within  10  to  100-fold  of  those  obtained  with 
conventional  sources,  and  linear  dynamic  ranges  were  about  four  decades. 
Progress  in  both  theory  and  experimental  achievements  has  been  rapid 
in  the  ensuing  years  (12,  33-38).  With  more  and  more  powerful  lasers 
available  (37,  39-45),  the  possibility  and  theory  of  non-linear 
phenomena  has  been  investigated  (40-44).  Several  authors  have  ex- 
perimentally achieved  near  or  complete  saturation  of  atomic  (46-50) 
and  molecular  systems  (51,  52). 


Laser  saturation  of  atomic  transitions  in  this  work  has  made 
possible  diagnostic  procedures  with  spatial,  temporal,  and  flame  com- 
position independence.  This  will  allow  analytical  chemists  to  measure 
the  quantum  efficiency,  total  number  density,  and  noise  power  density 
spectrum  in  an  atomizer  without  the  above  mentioned  problems  of  linear 
spectroscopy.  However,  the  benefits  of  this  freedom  must  be  weighed 
against  the  need  for  sophisticated  measurement  systems  (fast  time 
scale)  and  laser  temporal  and  spatial  inhomogeneity  effects. 


CHAPTER  2 
THEORETICAL  CONSIDERATIONS 

Quantum  Efficiency  for  a  Two-Level  Atom 

A  two-level  atom  is  one  which,  although  it  may  have  many  quantum 
states,  has  only  one  excited  state  available  for  thermal  population  in 
addition  to  the  ground  state.  After  it  has  been  delivered  to  the 
excited  state  by  absorption  of  a  photon,  it  can  be  returned  to  the 
ground  state  by  either  of  two  paths,  emission  of  a  photon  or  colli- 
sional  deactivation  (quenching).  Quantum  efficiency  (Y)  is  a  measure 
of  the  efficiency  of  the  radiative  process  with  respect  to  the  absorp- 
tion process.  A  schematic  diagram  of  a  three-level  atom  and  its 
processes  is  depicted  in  Figure  1.  A  two-level  atom  can  be  pictured 
by  ignoring  all  processes  connecting  the  third  level  and  the  other 
two  levels.  A  real  example  of  a  two-level  atom  is  snown  in  Figure  2. 
Only  the  4s  level  is  available  for  thermal  population  in  analytical 
flames  (temperature  range  from  about  2000  K  to  40C0  K,  k7  is  therefore 
0.2  to  0.37  eV). 

The  basic  fluorescence  radiance    expression  (53)  is  given  by: 


Bp=   (JL)  Y2lEVl27\dv  (1) 


where 

1  =  path  length  in  direction  of  detection  system,  m 

4ir  =  number  of  steradians  in  a  sphere  (fluorescence  is 
isotropic),  sr 


Figure  1.  Three  level  atom  model  photon  and  collisional 

processes  are  illustrated.  Symbols  are  explained 
below.  The  lower  state  is  symbolized  by  "1"  and 
the  upper  by  "u. " 

B,    Einstein  coefficient  of  induced  absorption 
lu   (n.3  j-1  s-I  Hz) 

B  i   Einstein  coefficient  of  induced  emission 
ul   (m3  j-1  s-1  Hz) 

A  ,   Einstein  coefficient  of  spontaneous  emission 
Ul   (s-1) 

Ev,   Spectral  irradiance  of  exciting  radiation 
lu  (W  m-2  Hz'l) 

k    Rate  constant  for  collisional  process  con- 
necting state  m  with  state  n 
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Figure  2.  Calcium  term  diagram  levels  1  and  2  are  indicated 
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Y«j  =  fluorescence  power  (quantum)  efficiency,  W  fluor- 
esced/W  absorbed 

L,j,9   =  spectral  irradiance  of  exciting  radiation  at  absorption 
iC       line,  Vl2,   W  m"2  Hz"1  (1  W  =  1  Js"1) 

J   k  dv  =  integrated  absorption  coefficient  over  absorption 
0      line,  m-1  Hz 

; 

c°° 

The  product  Ey1?  J   k  dv  is  the  power  absorbed  from  the  source  by  the 
analyte  atoms  per  meter  cubed  of  atomic  species.  The  integrated 
absorption  coefficient  is  given  by: 


hv1?        g.n-    , 


where 

hv,9  =  energy  of  the  exciting  photon,  J 

■1 


c  =  speed  of  light,  m  s 

B.o  =  Einstein  coefficient  of  induced  absorption,  m°  J 
u       s"1  Hz 


3  ,-1 


g..,  g?  =  statistical  weights  of  states  1  and  2,  respect"!  ve]y: 
dimension! ess 


n..,  n2  =  concentration  of  states  1  and  2,  respectively, 
{note  that  nj  +  n?  =  n-p  the  total  concentrate 


m 
ion  o 
atoms  in  all 'states' 


The  bracketed  quantity  corrects  for  the  effective  decrease  in  absorp- 
tion caused  by  stimulated  emission  from  the  upper  state. 

Using  the  rate  equation  approach,  and  invoking  steady  state,  we 
arrive  at 

(k12  ♦  B-^)   n,  ■  (kn   ♦  A21  +  bitty  ,z  (3! 

where 
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k,?,  k?,  =  excitation  and  deexcitation  non-radiational 
(collision)  rate  constants,  s~* 

Ap-,  =  Einstein  coefficient  of  spontaneous  emission,  s~ 

B0,  =  Einstein  coefficient  of  induced  emission,  m  J" 
21   s-1  Hz 

3 
B,9  =  Einstein  coefficient  of  induced  absorption,  m  J  ' 

u      s"1  Hz 

_3 

n, ,  n?  =  concentrations  of  electronic  states  1  and  2,  m 

c  =  speed  of  light,  m  s~ 


The  quantum  efficiency  is  defined  as 

A. 


21 


21   An  +  k21 


(4) 


and  A2i   is  related  to  EL-   and  B.p  by 

3  3 

8irhv  19  8irhv  10       g- 

A21  '  t-^2'  B21  '  (-T^)  *&   B12  W 

c  c      32 

Combining  these  expressions,  one  arrives  at 
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-2   -1 
where  E*vio  ""s  a  modified  saturation  spectral  irradiance  (W  m   Hz  ) 

evaluated  at  frequency  v.p  anc!  1S  defined  as 

CM21 

U        B21T21 


The  term  E*v.0  can  be  expressed  in  terms  of  the  saturation  spectral 
irradiance  E  v.p,  which  is  the  spectral  irradiance  required  to  bring 
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about  50%  of  the  maximum  fluorescence  radiance  oossible.  If  Es,,,n  is 


12 


expressed  as  a  function  of  E  v-2,  then 


s  9i 

F     =  F*    (        L 
v12     vi2  ^g1+  g2 


(8) 


Substituting  for  r^  in  terms  of  nT   (nT  =  n.   +  n?  and  using  Eqs.    3, 
6,  and  8) 
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The  maximum  fluorescence  is  then 
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When  EV12  =  E   2  then  Bp  =  Bp   /2, 
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A  theoretical   plot  of  Bp  vs.   Ev-2  appears   in  Figure  3.     If  a  plateau  in 
fluorescence  radiance  can  be  observed  as   laser  power  is  varied,  and   trie 
laser  power  at  50%  of  this  maximum  fluorescence  radiance  can  be  found, 
then  the  quantum  efficiency  can  be  found  from  the  quantity  ESV,?.      For 
calculations  in  this  work  the  quantity  ESx-2  will   be  used   (see  Eq. 
12). 
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The  necessary  conditions  for  the  measurement  of  quantum  efficiency 
in  this  way  are 

1.  Saturation  must  be  completely  achieved,  so  that  BF   can 

max 
be  evaluated. 

2.  One  must  measure  EVl2  by  measuring  the  following  quanti- 
ties: 

a.  laser  peak  power  (W) , 

2 

b.  laser  cross-sectional   area   (m  ), 

c.  laser  optical   bandwidth   (Hz). 

3.  Self-absorption  must  be  negligible. 

4.  The  source  must  be  a  spectral   continuum  with  respect  to 
the  absorption  line  width. 

5.  A  steady-state  condition  must  be  achieved. 

6.  No  coherence  effects  can  be  manifested.     The  rate  equation 
approach  must  be  valid. 

The  achievement  or  measurement  of  these  conditions  is  described  in 
Appendix  1. 

Quantum  Efficiency  for  a  Three-Level   Atom 

The  same  approach   (as  used  for  the  two-level   atom)  yields  for  a 
three-level   atom  like  Na   (see  Figure  4) 


V13      go 


1 


1  + 


'32 


v13 


(13) 


93       A23  +  k2l   +  k23 


where 


■18- 


'1/2 


1/2,   1   1/2 


11/2,2  1/2 


4s 


3s 


4p^ 


\  3d' 


i 


3d 


Figure  4.     Sodium  term  diagram  levels   1,   2,  and  3  are  indicated 
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k   =  collisional  rate  constant  connecting  atoms  in  the 
mn   m  state  with  those  in  the  n  state 

The  ESv-.3  is  again  related  to  Y-..,  but  no  simple  relationship 
exists.  The  collisional  constants,  which  are  unknown  in  the  litera- 
ture, interfere  with  the  interpretation. 

For  some  limiting  cases,  where  level  2  is  close  to  level  3  (e.g., 
sodium,  see  Figure  4)  or  is  close  to  level  1  (e.g.,  indium,  see 
Figure  5),  certain  assumptions  can  be  made. 

With  sodium,  for  instance,  if  the  assumption  is  made  that 


then 


k23  >:>  k21  +  A21 
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Then 
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E'vl3  '  <g  *  I]  *  ^ >  E*n3  <16' 


This  is  equivalent  to  a  two-level  system,  as  the  top  2  levels  have 
merged  into  one,  effectively. 
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Figure  5.  Indium  term  diagram  levels  1,  2,  and  3  are  indicated, 


-21- 

For  a  three-level  atom  like  In 
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The  same  type  of  assumption  should  hold  for  levels  1  and  2  as 
previously  for  level  2  and  3.  Therefore,  if 


then 


k21  >>  k32  +  A32 


k21       1      1 


and 


k21  +  k12  "  i  +^2  "    g2  "  q1   +  92 
K21      gl 
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E"vi3  =  (9l  +  g2  +  g3)  E*v13  (17] 


This  is  again  equivalent  to  a  two-level  system. 

Saturation 

Complete  saturation  means  that  the  population  of  the  excited 
state  and  the  ground  state  are  equalized  (with  due  regard  to  the  re- 
spective degeneracies).  Therefore,  the  kinetic  drains  on  the  excited 
state  are  negligible  with  respect  to  the  optical  pumping  rate.  A 
further  increase  in  irradiance  cannot  effect  an  increase  in  fluorescence 
radiance. 

A  measurement  of  this  (saturated)  state  usually  means  that  when 
the  irradiance  is  reduced,  via  some  suitable  filter,  the  fluorescence 


signal  does  not  decrease.  By  placing  calibrated  filters  in  the  beam, 
one  can  determine  the  point  at  which  only  half  of  the  signal  is  left; 
this  point  is  the  saturation  spectral  irradiance. 

Quantum  Efficiency  Via  the  Slope  Method 
Two-Level  System 

From  the  two-level   theory,   it  can  be  shown  that 
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A  plot  of  the  theoretical  shape  of  1/BF  vs.  l/E,  is  shown  in  Figure  6. 

The  pertinent  features  are  the  slope  and  the  intercept.  The  slope 
of  this  plot  contains  the  quantum  efficiency  (and  ru),  and  the  inter- 
cept contains  nT> 

For  this  procedure,  all  of  the  above  mentioned  criteria  must  be 
met  with  the  exception  of  completeness  of  saturation.  The  closer 
saturation  can  be  achieved,  however,  the  less  extrapolation  is 
necessary.  In  addition,  the  fluorescence  collection  depth,  1,  must 
be  measured,  and  the  collection  photometer  calibrated. 

If  both  sides  of  Eq.  18  are  multiplied  by  the  quantity 

4rr 
then  the  1/BF  axis  can  be  scaled  absolutely,  as  the  intercept  is 
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equal  to  the  quantity 


which  is  well  known. 


Three-Level  Case 


-25- 


'1  T  y2 


Sodium.     A  rearrangement  of  the  three-level   fluorescence  expres- 
sion yields 
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When  this  is  rearranged,  we  have 
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A  plot  of  these  quantities  cannot  be  scaled  absolutely  because  the 
intercept  includes  unknown  collisional  rate  constants.  The  same 
speculation  as  before  can  still  be  made,  however,  regarding  the  rate 
constants,  allowing  calibration  to  be  made  in  certain  limiting  cases, 
That  is,  k23  can  be  expected  to  be  much  larger  than  either  A01  or  k01. 
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The  term 
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then  reduces  to 
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By  the  law  of  mass  action  at  equilibrium 
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The  fluorescence  slope  method  equation  reduces  to 
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Rearranged  this  equation  resembles  that  for  two  levels 
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Indium.  For  an  atom  of  the  indium  type,  where  the  second  ex- 
cited state  is  very  close  to  the  ground  state, the  slope  method  equa- 
tion is 
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In  this  case  k21  can  be  assumed  to  be  much  larger  than  A-2  +  k~?: 
and  this  term  then  is  negligible. 
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Total  Number  Density 
Two-Level  Absolute  Plateau  Method 


A  measurement  of  Bf 

max 


h 


\ax  ■  te^,^^ 


produces  n-j.  easily  if  the  measurement  of  Bp    is  absolute  (see  Appendix 
1)  and  if  1,  A0       and  the  g's  are  known. 


max 


Two-Level-Slope  Method 


The  calibration  of  the  l/EL  axis  is  direct  for  total    number  density 
if  g^ ,   g.p ,  h,  \>.„t   Bp,   and  A21    are   known,   as  the  intercept  of  the 


1/B-  axis  is 


4ir  91  +  32 


B-   lhv12A2lnT    g2 

Three-Level  Case— Absolute  Plateau  Method 

Unless  assumptions  or  experiments  are  made   for  the  relative  size 
of  the  collisional   constants,  calibration  is  impossible.     Relative 
measurements   are  still   valid,  and  saturation  still   assures   immunity 
to  quantum  efficiency  effects  as  in  the  two-level   case. 

Three-Level   Case— Slope  Method 

Once  again  assumptions  must  be  made  for  collisional   constants   in 
order  to  permit  calibration  of  Bp  in  terms  of  ru. 

Possibility  of  Absolute  Calibration 

If  the  ratio  of  resonance  to  anti-stokes   fluorescence  is   taken, 
in  saturation,  we  have 
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Both  ratios 
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are  probably  not  much  greater  than  one,  as  k?.  is  expected  to  be  at 
least  comparable  to  either  numerator.  Then,  if  E«i  »  kT,  two  con- 
ditions are  satisfied 
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Equation  22  then  reduces  to 
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If  a  flame  of  known  temperature  can  be  produced,  then  the  ratio 
^A32+k32^k21  can  ^e  evaluated.  This  allows  Y  (and  n-J  to  be  obtaine 
in  absolute  units.     For  a  flame  with  a  temperature  of  2300  K  and 
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thallium  as  a  probe  (E.2  =  0.966  eV),the  exponential  term  is  as 
follows: 


exp[+0. 966/2300  x  0.861  x  10"4]  =  131 


For  thallium 


Then 
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!F^3     2  +  A32  +  k32 
2+3max  .       k2l 


BF  135 

l-e-3 

:  ,max 
l-»3 


The  difference  between  the  ratio  and  2/135  =  0.0148  is  the  contribu- 
tion of  A-2  +  kop/koi*  This  value,  when  inserted  into  Eq.  21  (for 
measurement  of  thallium,  in  this  case),  would  allow  the  intercept  of 
the 

nT1A31hvl3 
4ttBp 
rl«-3- 
1+3 

axis  to  be  completely  calibrated,  as  the  intercept  then  contains  only 
known  quantities.  If  BF  is  known  absolutely  (calibrated),  then  nT  can 
also  be  arrived  at  unambiguously. 

Noise  Power  Density 

The  theory  and  practical  measurement  of  power  spectra  (power  vs. 
frequency  data)  have  been  adequately  covered  by  Blackman  and  Tukey  (54). 
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Two  approaches  are  possible  to  measure  noise  power  spectra.  The  first 
is  to  compute  the  autocovariance  function  for  the  data  collected.  The 
Fourier  transform  of  this  is  the  power  spectrum.  The  second  procedure 
is  to  compute  the  Fourier  transform  of  the  data,  then  to  square  and 
add  the  frequency  and  phase  components  to  find  the  power  spectrum. 
Although  both  methods  were  available,  the  second  one  was  chosen  for 
its  conceptual  familiarity.  The  use  of  a  discrete  data  collection 
system  (analog  to  digital  converter)  required  a  discrete  Fourier 
transform,  with  attendant  aliasing  and  bandwidth  problems. 

The  application  of  saturated  atomic  fluorescence  will  result  in 
an  ability  to  separate  the  noises  resulting  from  the  nebulization  pro- 
cess from  those  stemming  from  the  atomization  of  the  analyte.  Since 
these  two  sources  of  noise  are  independent,  their  noise  power  spectra 
add  linearly.  Since  the  measure  being  computed  in  our  case  is  the 
noise  current  spectrum,  these  would  add  quadratically.  If  an  analyte 
such  as  copper,  which  is  completely  atomized  in  a  standard  flame,  is 
measured,  the  noise  power  under  known  conditions  can  be  found. 

This  noise  power  should  be  reproducible  in  nature  with  the  ex- 
ception of  certain  non-stationary  noise  sources.  If  this  proves  to  be 
the  case  then  the  total  noise  and  the  1/f  noise  pertinent  to  most 
spectroscopic  measurements  can  be  investigated  with  respect  to  design 
elements  and  operating  characteristics  common  to  atomic  spectroscopy. 
These  include  burner,  chamber,  torch,  and  nebulizer  design,  as  well  as 
gas  flows  and  flame  stoichiometry.  Noises  involved  in  the  flame  or 
plasma  processes  up  to  the  point  of  atomization  could  then  be  sub- 
tracted from  the  noises  seen  for  an  incompletely  atomized  analyte, 
producing  the  noise  inherent  in  the  atomization  process. 
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The  unfortunate  truth,  however,  is  that  many  of  the  underlying 
noise  sources  in  analytical  spectroscopy  are  non-stationary. 
Examples  of  non-stationary  noises  are  nebulizer  clogging,  RF  noise, 
60  Hz  mains  noise,  some  source  drift  noises,  stray  light  noise  (room 
and  daylight),  and  some  electronic  drift  noise.  These  are  not  really 
approachable  via  the  technique  or  theory  of  noise  power  spectra,  though 
often  a  reasonable  idea  of  the  noises  can  be  found  (55). 


CHAPTER  3 
EXPERIMENTAL 

General  Comments 

A  variety  of  experimental  systems  were  used  for  collection  of  data 
pertinent  to  these  measurements.  The  basic  system  shown  in  Figure  7 
is  the  laser  excited  atomic  fluorescence  system  described  by  Weeks 
et  al .  (13).  The  system  has  been  modified  slightly,  as  reported  by 
Bower  et  al .  (56)  and  Bradshaw  et  al .  (10).  Several  of  the  features 
changed  since  the  measurements  of  Weeks  et  al.  are  of  high  importance. 
The  spatial  effects  of  saturation  in  the  wings  of  focussed  Gaussian 
laser  beams,  as  observed  by  Blackburn  et  al.  (50),  and  previously 
predicted  by  Rodrigo  (57)  and  theoretically  approached  by  Daily  (58) 
have  been  eliminated  by  moving  the  observation  point  approximately  10 
feet  further  from  the  laser.  At  this  point,  the  spatially  homogeneous 
section  of  the  laser  beam  was  enveloping  the  whole  seeded  flame  (see 
Appendix  1). 

In  order  to  facilitate  alignment,  it  was  necessary  to  add  folding 
mirrors  (M)  with  adjustments  in  the  horizontal  and  vertical  planes. 
A  problem  involving  dye  changes  leading  to  a  different  angle  of 
emergence  has  been  circumvented    by  the  use  of  the  folding  mirrors 
and  the  combined  aperture-scatter  shields  (ASS).  The  two  apertures 
initially  were  set  to  define  the  beam  line  of  interaction  of  the  dye 
laser  and  the  flame.  When  subsequent  dye  changes  produced  non-col1near 
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Table  1 
Laser  Excited  Atomic  Fluorescence  Equipment  List 


Item 


Manufacturer 


Model  UV-14  Nitrogen  Laser 
Model  DL-300  Tunable  Dye  Laser 
Model  FL-2000  Tunable  Dye  Laser 
Nebul izer  chamber 
Capillary  burner 
Flame-shielded  burner 
Model  H-10  monochromator 

Model  412B  High  Voltage 

Power  Supply 

Model  R106  Photomultiplier 
Tube 

Model  160  boxcar 

Model  162  boxcar  mainframe 

Model  163  boxcar  plug-in 
module 

Model  164  boxcar  plug-in 
modul e 

Type  S-2  sampling  heat 

Servoriter  strip  chart  recorder 

Apertures  (iris  diaphragms) 


Molectron  Corp.,  Sunnyvale,  Cal. 

Molectron  Corp.,  Sunnyvale,  Cal. 

Lambda  Physik,  Gottingen,  W.  Ger. 

Perkin-Elmer  Corp.,  Norwalk,  Conn. 

laboratory  built 

laboratory  built 

Jobin-Yvon,  %   ISA,  Metuchen,  N.J. 

John  Fluke,  Seattle,  Wash. 

Hamamatsu  Corp.,  Middlesex,  N.J. 

Princeton  Applied  Research, 
Princeton,  N.J. 

Princeton  Applied  Research, 
Princeton,  N.J. 

Princeton  Applied  Research, 
Princeton,  N.J. 

Princeton  Applied  Research, 
Princeton,  N.J. 

Tektronix,  Portland,  Oregon 

Texas  Instruments,  Dallas,  Texas 

Edmund  Scientific,  Barrington,  N.J. 
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beams,  these  could  be  adjusted  with  the  mirrors  to  pass  through  the 
apertures. 

The  collection  system  was  aligned  at  a  90  degree  angle  to  the  in- 
coming beam  by  the  use  of  helium-neon  alignment  laser  crossed  with  the 
dye  laser  beam. 

Radio  Frequency  Shielding 

Radio  Frequency  (RF)  shielding  of  the  entire  Molectron  nitrogen 
pump  laser  enabled  the  output  to  be  freed  of  the  large  drift-type  noise 
associated  with  people  and  objects  moving  in  proximity  to  the  experi- 
mental system.  This  involved  enclosing  the  whole  pump  laser  and  power 
supply  in  a  brass  screen  Faraday  case.  All  electronics  were  connected 
via  grounding  straps  to  a  salt-bed  earth  ground.  All  cables  leading 
from  the  photomultipl ier  to  the  detection  electronics  were  double- 
shielded  by  a  coaxial  woven  shield  augmenting  the  internal  shield. 
These  modifications  enabled  RF  noise-free  measurements  into  the  milli- 
volt level,  which  corresponded  (with  a  50  SI  2-fold  attenuator)  to 
40  yA  of  peak  current.  This  was  possible  even  with  a  background 
current  of  100  yA.  For  most  measurements,  the  limiting  noise  was  the 
electronic  noise  of  the  boxcar  averaaer. 
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Photomu 1 ti  pi i  er 

The  photomultiplier  was  wired  for  fast  pulse,  high  current 
output,  as  reported  by  Fraser  and  Winefordner  (32)  and  operated  at 
-1000  V.  Cable  of  50  n  impedance  (RG-58U)  was  used  for  all 
connections. 


Fluorescence  Flux  Collection 

A  Jobin-Yvon  H-10  low  dispersion  monochromator  was  used 
for  collection  of  fluorescence.  Dispersion  is  low  and  collection 
efficiency  high,  as  suggested  by  Weeks  et  al.  (13).  A  two  Tense 
collection  system  with  an  intervening  aperture  was  used  to  help 
cut  down  on  scatter  for  resonance  collection  cases.  Scatter  in 
these  cases  was  reduced  by  liberal  use  of  light  traps  (LT),  aper- 
ture scatter-shields  (ASS),  and  black  felt  cloth.  In  many  cases, 
this  scatter  signal  was  Rayleigh  scatter  (indicated  by  an  increase 
when  the  flame  was  turned  off),  and  could  not  be  completely  elimi- 
nated. In  such  cases,  the  noise  limit  for  measurements  was  due  to 
laser  carried  (peak  to  peak  variation)  scatter  noise.  The  image 
formed  by  the  two  lenses  at  the  slit  of  the  monochromator  was  of 
unit  magnification. 
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Detection  Electronics 

A  PAR  162  boxcar  averager  was  used  for  detection  and  measurement 
of  fast  fluorescence  pulses.  Two  different  inserts  (plug-ins)  were 
used  in  measurement  and  alignment.  The  164  plug-in  offered  fixed 
sensitivity  (100  mV  full  scale)  and  several  integration  times.  It  was 
quite  hardy,  and  was  used  for  preliminary  alignment,  when  a  very  large 
scatter  signal  was  used  to  set  the  gate  delay.  The  50  a   input  im- 
pedance was  always  used  to  match  the  cable  and  minimize  ringing. 
Some  ringing  was  still  apparent,  however.  The  163  plug-in  with  an  S-2 
sampling  head  was  used  for  the  main  measurement  scheme.  The  gate  was 
75  ps  (S-2  sampling  head);  jitter  in  the  boxcar  gating  circuits  led  to 
some  distortion  of  this  aperture.  This  head  was  used  to  find  plateau 
(saturation)  values  on  top  of  the  fluorescence  temporal  pulses,  which 
corresponded  to  saturation  of  the  atomic  fluorescence. 

The  synchronous  trigger  of  the  laser  could  not  be  used,  as  RF 
feedback  occurred  at  the  time  of  each  trigger  and  interfered  with  data 
collection.  A  separate  trigger  generator  was  used,  which  had  a  delayed 
pulse  feature.  This  enabled  correct  timing  between  laser  firing  and 
boxcar  gating. 

Nitrogen  Laser 

The  nitrogen  pump  laser  was  run  according  to  manufacturer's  in- 
structions. The  high-voltage  power  supply  was  kept  at  22  kV  to  give 
better  longevity  of  the  critical  components  (thyratron  and  capacitor), 
as  these  had  been  problematical  in  the  past.  The  nitrogen  flow  rate 
was  15  liters  per  minute  for  all  experiments.  Operating  pressure  was 
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maintained  at  55  torr.  The  pressure  transducer  had  to  be  replaced 
and  calibrated  twice,  however.  This  involved  approximately  half  a 
day  to  a  day's  work,  since  the  entire  pump  head  had  to  be  removed  from 
the  Faraday  cage.  All  connections  had  to  be  made  again  outside  the 
cage  in  order  to  test  fire  and  calibrate  the  unit. 

The  repetition  rate  was  kept  at  about  15-20  Hz  for  all  experi- 
ments. Greater  repetition  rates  led  to  lower  peak  powers.  High  peak 
power  was  necessary  for  saturation.  However,  a  rate  below  ca.  20  Hz 
would  lead  to  signal  leakage  problems  on  the  boxcar  mainframe. 

An  interesting  phenomenon  was  observed  at  the  20  Hz  repetition 
rate.  A  beat  frequency  appeared  on  the  dye  laser  output.  This  was 
first  attributed  to  the  power  supply  of  the  laser,  as  its  voltage 
dipped  when  the  laser  was  fired.  This  idea  was  discarded  when  it  was 
noticed  that  the  beat  frequency  was  about  2  Hz.  The  observation  was 
made  that  as  the  dye  fluoresced  in  the  cuvette  on  each  nitrogen  laser 
pulse,  the  magnetic  stirrer  appeared  to  be  precessing  at  about  2  Hz. 
The  motor  provided  was  an  AC  motor  running  at  a  multiple  of  60  Hz. 
The  stirring  of  the  dye  cell  seemed  to  be  incomplete  enough  that  when 
the  nitrogen  pulse  came  at  a  sub-multiple  of  the  stirring  rate,  a 
beat  frequency  due  to  excitation  of  spent  dye  appeared.  A  need  for 
improvement  in  stirring  was  quite  apparent.  This  beat  frequency 
effect  had  never  appeared  in  the  literature. 

Voltage  drifts  during  the  day  were  so  bad  (1  V  on  regulated 
supplies,  10  V  on  unregulated  lines)  that  severe  trigger  drift  became 
apparent.  The  nitrogen  pump  laser  was  yery   susceptible  to  trigger 
drifts  through  its  high  voltage  threshold  performance.  A  voltage 
change  of  1  V  on  the  120  VAC  line  can  induce  a  20  ns  drift  in  the 
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trigger  to  firing  delay.  Using  a  164  plug-in  on  the  boxcar,  this  was 
annoying.  When  a  163  plug-in  is  used  the  signal  can  entirely  dis- 
appear inside  of  5  min  time.  Therefore,  all  of  the  experiments  had 
to  be  done  at  night  when  stability  was  much  better,  but  still  trouble- 
some, however.  Generally  during  each  series  of  measurements  the 
trigger  had  to  be  reset  at  least  once. 

An  optical  trigger  was  used  for  subsequent  experiments.  The 
163  plug-in,  however,  was  not  used.  Indications  were  that  this 
arrangement  was  much  more  stable  with  respect  to  jitter,  especially 
if  triggered  from  the  pump  laser.  Use  of  the  optical  trigger  necessi- 
tated placing  a  delay  line  (V155Z050,  Allen  Avionics,  Mineola,  NY) 
into  the  signal  line  from  the  photomultiplier  to  produce  a  75  ns  delay 
necessary  to  compensate  for  the  internal  gate  delay  of  the  boxcar. 
A  50  q   attenuator  was  placed  before  the  delay  line  to  produce  a 
voltage  pulse  from  the  current  pulse.  Indications  were  that  distor- 
tion and  noise  introduction  by  the  delay  line  were  minimal. 

Dye  Laser  Operation 

The  dye  laser  was  also  operated  according  to  manufacturers  in- 
structions. Dyes  were  prepared  according  to  Table  2.  Adjustments  of 
the  dye  cell  carriage  were  made  after  each  new  dye  cell  was  introduced 
into  the  beam  line.  This  should  not  have  been  necessary  if  the  laser 
was  functioning  correctly.  There  was  some  indication  that  this  was 
an  artifact  due  to  the  state  of  our  dye  cells.  Possibly  the  anti- 
reflective  coating  had  worn  off  in  several  years  of  use  (see  Appen- 
dix 1). 
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Table 

2 

Laser 

Dye  Preparation 

Dye 

Concentration 
(M) 

Solvent 

Wavelength 
Range  (nm) 
(10%  points)* 

DPS 

Saturated 
1.2  x  10"3 

p-dioxane 

396-416 

Bis-MSB 

1.2  x  10"3 

p-dioxane 

411-430 

7D4MC 

ID"2 

ethanol 

440-478 

R6G 

5  x  10"3 

etnanol 

568-605 

*The  10%  points  are  the  wavelengths  at  which  laser  output  is  10% 
of  the  peak  output. 
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Dye  was  changed  frequently  to  avoid  power  loss  due  to  dye  deteri- 
oration. This  was  done  in  one  of  two  ways.  The  method  of  exchange 
recommended  by  the  manufacturer  was  to  remove  the  dye  via  a  pi  pet, 
rinse  the  cell  with  fresh  dye,  and  refill  with  fresh  dye.  This  was 
moderately  successful.  Some  dyes  (7D4MC  in  particular)  seemed  to 
function  better  when  the  dye  cell  was  removed  entirely  from  the  dye 
head,  rinsed  thoroughly  with  alcohol,  then  fresh  dye,  and  then  re- 
filled with  fresh  dye.  Peak  power  with  7D4MC  exchanged  in  this  manner 
was  enhanced  at  least  2-fold  over  the  more  cursory  method  of  exchange. 
Bubbling  of  pure  nitrogen  for  two  minutes  through  the  dyes  (as  sug- 
gested by  Exciton)  did  not  seem  to  make  a  substantial  change  in  dye 
longevity.  Perhaps  nitrogen  saturation  during  preparation  and  storage 
would  have  helped  more. 

Flame  System 

A  gas-flow  stabilization  system  consisting  of  individual  pressure 
gauges  and  rotameters  was  used  for  each  gas  flow.  The  entire  system 
was  calibrated  by  means  of  a  linear  mass  flow  meter  (ALK-50K,  Hastings, 
Hampton,  Virginia).  Each  time  a  flow  was  measured,  the  pressure  was 
readjusted,  and  the  flow  and  rotameter  readings  recorded.  The  com- 
bination of  a  fixed  pressure  and  the  rotameter  readings  led  to  a  very 
reproducible  and  precise  system.  Premixing  of  flame  gases  was  done 
to  avoid  background  noise  due  to  poor  mixing  in  the  chamber  of  the 
burner.  The  flows  passed  into  the  appropriate  ports  of  a  nebulizer- 
flow  chamber  (Perkin-Elmer  Corporation,  Norwalk,  Conn.).  The  same 
composition  (premixed)  passed  into  a  separated  shielding  flame  for 


hydrogen-based  flames.  The  fog  and  premixed  combustion  gases  left 
the  flow  chamber  and  passed  into  a  capillary,  flame-shielded  burner 
for  hydrogen  based  flames  (previously  described  by  Snelleman  (59)). 
There  they  were  combusted.  The  flame-shielded  flame  was  surrounded 
by  an  inert  gas  sheath  to  avoid  background  DC  emission  and  noise  as 
much  as  possible. 

Acetylene  based  flames  were  combusted  above  a  capillary  burner  (Haragu- 
chi  and  Winefordner  (60))  supported  on  the  same  nebulizer  chamber. 
Flows  for  the  different  flames  used  are  shown  in  Table  3. 

Solutions 

Stock  solutions  of  1000  ug/ml  for  all  elements  were  made  from 
reagent  grade  chemicals  in  deionized  water  per  Parsons  et  al .  (61) 
(see  Table  4). 

Measurement  Procedure  for  Y  and  nT 

The  necessary  criteria  for  evaluation  of  quantum  efficiency  and 
total  number  density  via  the  two  different  schemes  are 

1.  Saturation  must  be  complete  (for  plateau  method)  or  nearly 
complete  (for  the  slope  method);  that  is,  an  increase  in  laser 
power  must  not  result  in  an  increase  in  fluorescence 
radiance. 

2.  One  must  measure  laser  spectral  irradiance  via  measurement 
of  the  following  quantities: 

a.  laser  optical  bandwidth--via  some  type  of  wavelength 
scan, 
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Table  3 
Flame  Gas  Flow  Rates 


Air-Acetylene  Flame 

Air         9.88  1/min 

Acetylene    1.65  1/min 

Hydrogen-Oxygen- Inert  Gas  Flame 

Hydrogen  1.62  1/min 

Oxygen  0.81  1/min 

Nitrogen  or  Argon  4.51  1/min 


-45- 


Table  4 
Source  of  Reagents 


Sodium  NaCl  Mallinckrodt  Chemical  Works, 

St.  Louis,  Mo. 

Calcium  CaCCL  Mallinckrodt  Chemical  Works, 

St.  Louis,  Mo. 

Strontium  SrC03  Fisher  Scientific, 

Fair  Lawn,  N.J. 

Indium  Ir^CL  Apache  Chemicals,  Seward,  111 
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b.  laser  peak  power— via  calibrated  photosensitive  device, 

c.  laser  cross-sectional  area—via  geometrical  considera- 
tions; b  and  c  measurements  may  be  combined. 

3.  Self-absorption  must  be  negligible.  (This  is  automatically 
achieved  if  saturation  is  complete,  as  the  absorption  coef- 
ficient then  goes  to  zero.) 

4.  The  source  must  be  a  (pseudo-)continuum  across  the  line- 
width  of  absorption;  i.e.,  the  wavelength  spread  of  the  laser 
must  be  measured  and  determined  not  to  vary  substantially 
across  the  absorption  linewidth. 

5.  Steady  state  fluorescence  must  be  achieved.  The  temporal 
behavior  of  the  atomic  fluorescence  pulse  must  be  observed. 
Measurement  of  the  saturated  value  must  be  during  a  suitable 
steady  state. 

6.  The  rate  equation  approach  must  be  valid.  No  coherence 
effects  may  be  seen. 

7.  The  fluorescence  depth  1  must  be  found  (Figure  8). 

In  addition,  if  the  possibility  of  absolute  values  of  Y  and  nT 
is_  to  be  evaluated  for  a  three-level  atom,  the  ratio  of  resonance  to 
stokes  fluorescence  must  be  found  in  a  flame  with  a  known  temperature. 

Since  the  proof  of  validity  of  these  criteria  is  somewhat 
tedious,  it  will  be  assigned  to  Appendix  1,  with  the  exception  of  the 
first  point,  completeness  of  saturation,  because  this  is  the  basis 
of  the  whole  measurement  scheme. 

For  each  dye  change,  alignment  of  the  laser  beam  with  the  burner 
head  was  undertaken.  Side  to  side  position  of  the  laser  was  checked 
with  the  flame  off,  by  the  use  of  the  aperture-scatter  shields  and  a 
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card  over  the  burner.  Vertical  positioning  was  accomplished  by  using 
the  card  to  project  an  image  of  the  laser  beam  through  the  first  col- 
lection lens  onto  the  slit-shaped  aperture.  After  spatial  alignment 
was  completed,  the  gate  delay  was  peaked  up  on  a  scatter  signal. 
Next,  the  flame  was  lit,  and  analyte  was  introduced.  The  dye  laser 
was  wavelength  scanned,  using  the  scan  control  until  a  substantial 
signal  was  found.  Since  saturated  atoms  show  a  broadened  response 
to  wavelength  (see  Appendix  1),  a  calibrated  neutral  density  filter 
was  set  in  an  accessory  optics  holder  in  the  laser  beam.  In  this  way, 
a  linear  response  to  laser  illumination  was  ensured.  Neutral  density 
filters  of  about  2.0  to  3.0  were  generally  required.  When  illumina- 
tion was  in  the  linear  region,  the  dye  laser  output  could  be  easily 
centered  on  the  wavelength  of  the  atomic  transition.  When  the 
neutral  density  filter  was  removed,  the  atoms  were  saturated. 
Generally  analyte  concentration  was  chosen  to  yield  a  high  enough 
signal  so  that  a  linear  portion  of  the  laser  fluorescence  vs.  power 
plot  could  be  reached  as  the  laser  power  was  decreased.  In  some 
cases,  this  involved  using  a  higher  concentration  at  low  laser  power 
to  ensure  sufficient  signal  to  noise  ratio.  Linearity  of  fluorescence 
with  respect  to  concentration  was  always  checked  when  this  was  done. 
The  upper  limit  of  concentration  was  set  to  avoid  problems  of  vapor- 
ization (especially  prevalent  in  the  low  temperature  hydrogen -based 
flames  used  here)  and  self-absorption.  No  absorption  was  noticed, 
however,  with  solutions  up  to  1000  yg/ml  in  concentration,  although 
the  measurement  of  absorption  was  complicated  by  severe  shot  to  shot 
peak  power  variations  on  the  dye  laser.  In  most  cases,  concentra- 
tions were  in  the  1  to  20  yg/ml  range. 
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The  procedure  for  measurements  started  with  a  check  for  PMT 
linearity.  Even  with  the  use  of  a  low  resistance  (high  current) 
dynode  chain  and  capacitor  charge  storage  for  each  dynode,  the  PMT 
output  current  could  become  non-linear  with  respect  to  illumination 
when  the  light  levels  became  high  enough.  Each  time  data  were  taken, 
a  0.3  neutral  density  filter  was  inserted  in  front  of  the  mono- 
chromator.  If  the  signal  did  not  decrease  by  50%,  a  lower  concen- 
tration was  aspirated,  and  checked  for  linearity.  Approximately  half 
of  the  full  scale  (1  V)  was  usable  for  the  163-S-2  plug-in  combina- 
tion. This  corresponded  (with  a  50  a   2-fold  attenuator)  to  a  10  mA 
peak  current.  This  compared  with  a  limit  of  0.3  mA  peak  current  limit 
of  Olivares  (46),  for  a  220  kft  resistor  chain  with  no  added  capaci- 
tance, for  a  1P28  PMT.  A  linear  dynamic  range  of  30-fold  higher 
(based  on  current)  was  obtained  in  our  case. 

After  the  linearity  check,  the  actual  saturation  curve  measure- 
ment was  performed.  This  consisted  of  measuring  full  scale  fluorescence, 
blank  scatter,  and  several  reduced  fluorescence  and  scatter  values 
when  calibrated  neutral  density  filters  were  placed  in  the  laser  beam 
singly  or  in  combination.  Scatter  corrections  for  resonance  measure- 
ments were  generally  significant  until  the  laser  intensity  was  re- 
duced 1  to  2  orders  of  magnitude  by  neutral  density  filters. 

These  fluorescence  measurements  were  converted  to  absolute 
radiance  values  following  calibration  of  the  optical  collection 
system,  the  boxcar,  and  the  recorder  used  to  collect  data  (see 
Appendix  1). 

From  the  absolute  radiance  measurements,  values  were  computed 
for  quantum  efficiency  and  total  number  density.  In  some  cases  rela- 
tive data  only  were  computed  for  different  flames. 
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Noise  Power  Density 

The  measurement  of  noise  power  density  spectra  in  academic  analy- 
tical chemistry  is  generally  accomplished  via  the  discrete  fast 
Fourier  transform  (DFFT).  This  approach  is  necessitated  by  the  lack 
of  continuous  measurement  power  spectrum  analysis  instrumentation  in 
most  laboratories.  The  dedicated  instrumentation  required  by  this 
method  is  generally  too  inflexible  to  be  affordable  under  tight 
budgeting  restrictions  of  this  time.  On  the  other  hand,  the  general 
availability  of  digital  computers  and  analog  to  digital  data  collec- 
tion systems  makes  discrete  (or  sampled)  data  collection  fit  reason- 
ably within  the  price  and  flexibility  range  of  most  laboratories. 
This  leaves  one  with  only  the  DFFT  approach  to  noise  power  spectra 
evaluation  possible.  This  enforced  choice  leads  to  a  loss  of  flexi- 
bility in  one  way,  and  to  a  gain  in  another  direction.  The  disadvan- 
tages of  discretely  sampled  records  for  approximating  a  power  spectrum 
are  very  well  covered  in  Blackman  and  Tukey  (54)  and  include  aliasing 
(necessitating  very  well  designed  filters  for  spectra  with  wide  dynamic 
ranges  and  wide  frequency  ranges),  intermodulation  distortion,  and 
finite  bandwidth. 

The  flexibility  advantages  are   those  which  stem  from  the  com- 
puter itself: 

1.  ease  of  collection, 

2.  ease  of  calibration, 

3.  fast  data  reduction, 

4.  ease  of  presentation, 

5.  flexibility  of  handling  of  data  and  power  spectra. 
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The  primary  purpose  of  noise  power  spectra  collection  in  this  work 
was  to  establish  a  background  for  later  possible  saturated  noise  power 
spectra,  which  should  show  the  advantages  enumerated  in  the  Introduc- 
tion section. 

The  primary  sources  of  noise  for  atomic  spectroscopy  consist  of 
the  following: 

1.  Source  related  noise 

a .  wh  i  te , 

b.  whistle  or  proportional, 

c.  flicker  or  1/f. 

2.  Background  noise-non-source  related 

a .  wh  i  te , 

b.  whistle  or  proportional, 

c.  fl icker  or  1/f. 

3.  Sample  related  noise 

a .  wh  i  te , 

b.  whistle  or   proportional, 

c.  flicker  or  1/f. 

The  preliminary  experiments  presented  here  measured  the  latter  two 
types  of  noises.  The  first  type  will  be  the  subject  of  later  experi- 
ments. 

Background  noise  can  be  somewhat  artificially  separated  into  two 
types.  The  first  type  is  related  to  the  bulk  properties  of  the  flame 
gases;  this  type  could  be  exemplified  by  the  emission  from  the  C? 
molecule,  a  product  of  partially  completed  combustion.  This  molecule 
has  a  major  emission  at  516.5  nm  called  the  Swan  band.  While  it  does 
not  interfere  with  most  atomic  transitions,  other  bands  can  interfere 
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with  DC  and  fluorescence  measurements  with  pulsed  sources  .  The 
Swan  band  is  sensitive  and  can  be  used  as  an  indicator  of  formation  of  C«, 
and  therefore  as  an  indicator  of  combustion  completeness.  This 
molecule  has  been  used  as  a  measure  of  success  for  premixing  com- 
bustion gases  (29).  Noise  power  spectra  collected  under  fixed  con- 
ditions also  provide  a  useful  figure  of  merit  for  the  evaluation  of 
success  of  mixing. 

A  second  type  of  background  emission  noise  source  is  one  which  is 
more  indicative  of  external  conditions.  An  example  is  the  OH  molecule, 
which  is  strongly  indicative  of  the  success  of  sheathing  a  flame. 
The  pertinent  (and  sensitive)  band  is  centered  on  306  nm,  and  extends 
to  the  red  far  enough  to  interfere  severely  with  copper  and  silver 
measurements   while   a   companion  band  at  285  nm  interferes  with 
lead  and  magnesium  measurements.  Noise  power  spectra  at  this  band 
can  hopefully  lead  to  better  design  and  construction  for  flame 
sheathing  systems. 

Sample  related  noise  can  be  evaluated  as  a  preliminary  to  satura- 
tion measurements  by  using  atomic  emission.  In  order  to  collect  data 
on  both  nebulization  and  atomization  noises  without  source  noises 
interfering,  it  is  necessary  to  measure  the  noise  involved  in  atomic 
emission,  unless  saturated  measurements  lead  to  an  improvement  in 
signal  to  noise  ratio.  An  atom  must  be  chosen  which  has  strong  emis- 
sion, and  is  not  near  a  background  spectral  region  which  would 
interfere.  Strontium  is  a  likely  candidate.  Its  emission  is  in  a 
relatively  background-free  region,  and  should  therefore  yield  infor- 
mation on  nebulization  and  atomization  efficiency  noise.  Strontium 
is  quite  sensitive  to  flame  temperature,  and  therefore  composition. 
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It  should  be  a  good  candidate  for  the  proposed  separation  of  information 
from  the  two  sources  of  noise.  In  addition,  it  is  among  the  most 
strongly  saturated  atoms  found  in  this  work.  The  general  scheme  of 
investigation  is  as  follows.  The  nebulization  chambers  and  nebulizers 
of  three  companies  are  to  be  investigated  within  the  guidelines 
established  above.  Conditions  to  be  varied  are  use  of  an  inert  gas 
sheath,  premix  chamber,  and  removal  of   the  nebulizer.  In  addition, 
data  collection  will  be  made  in  three  frequency  ranges  (0-50  Hz, 
0-500  Hz,  0-10,000  Hz). 


CHAPTER  4 
DATA 

Notation 

A  general  notation  system  is  developed  here  for  the  reporting  of 
the  experimental  results  in  this  chapter.  Ratios  will  be  listed  as  "R," 
with  addition  of  one  or  more  of  the  following  subscripts  for  clarifica- 
tion of  the  meaning: 

HON   Hydrogen-oxygen-nitrogen 

HOA   Hydrogen-oxygen-argon 

S    Slope  method  of  calculating  a  quantum  efficiency  or  ratio 
of  quantum  efficiencies 

E    E..  method  of  calculating  a  quantum  efficiency  or  ratio  of 
quantum  efficiencies 

B    BFmax  me'tnoci  °^   calculating  a  total  number  density  or  ratio 
of  total  number  densities 

I    Intercept  method  of  calculating  a  total  number  density  or 
ratio  of  total  number  densities 

A    Measurements  made  in  the  argon  diluted  hydrogen  based  flame 

N    Measurements  made  in  the  nitrogen  diluted  hydrogen  based 
f  1  ame 

Y    Quantum  efficiency 

n-j.   Total  number  density 

References  to  the  use  of  these  methods  can  be  found  in  Chapter  2 
(Theoretical  Considerations).  Ratios  "R"  will  always  refer  to  argon 
diluted  flame  values  over  nitrogen  diluted  flame  values. 
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Strontium 
Boxcar  Survey 

The  need  to  use  the  most  versatile  and  least  complicated  equipment 
possible  was  recognized  at  the  commencement  of  this  work.  One  of  the 
goals  of  the  work  was  to  find  the  optimal  instrumentation  for  flame 
diagnostic  measurement  via  saturated  atomic  fluorescence,  choosing  from 
the  available  choices  in  this  laboratory.  A  series  of  measurements  were 
made  under  set  conditions  for  the  comparison  of  three  boxcar  averager 
measurement  systems.  The  first  pair  of  measurements  (Figure  9  )  shows 
saturation  curves  under  identical  illumination  conditions  for  the 
Princeton  Applied  Research  (PAR)  160  boxcar  and  the  PAR  162  mainframe 
with  the  164  plug-in.  The  two  curves  are  very  nearly  identical,  though 
the  scatter  on  this  data  was  fairly  great.  (A  substandard  burner  system 
was  used  for  these  two  preliminary  curves.  It  was  replaced  after  these 
two  experiments  with  the  system  described  in  the  text.)  The  gates  on 
both  boxcars  were  the  same  (15  ns).  The  10  ps  stretch  feature  of  the 
160  boxcar  did  not  seem  to  affect  the  extent  of  saturation.  The  second 
pair  of  measurements  (Figure  10)  shows  a  comparison  of  the  160  boxcar 
with  the  162  mainframe  using  a  163  plug-in  and  a  S-2  sampling  head  (75  ps 
aperture).  In  addition  to  the  boxcar  change,  a  change  of  dye  was  per- 
formed. As  expected,  the  extent  of  saturation  improved.  The  shorter 
gate  "viewed"  mainly  the  temporal  plateau  section  of  the  fluorescence 
pulse  (see  Appendix  1),  while  the  longer  gate  averaged  substantially 
unsaturated  portions  of  the  pulse  with  the  saturated  portion.  The  higher 
power  output  of  the  new  dye  probably  accounted  for  some  portion  of  the 
behavior. 
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Based  on  the  results  of  these  measurements  the  rest  of  the  measure- 
ments were  made  with  the  162  mainframe  163  plug-in  combination  to  achieve 
an  accurate  representation  of  the  saturation  condition. 

Flame  vs.  Inductively  Coupled  Plasma  (ICP) 

An  acetylene-air  flame  and  an  ICP  were  compared  to  determine  whether 
saturation  resulted  in  information  about  quantum  efficiency.  Two 
parameters  were  evaluated  under  equal  illumination  in  the  flame  and  ICP, 
namely  the  saturation  curve  and  the  saturation  broadening  (Appendix  1) 
for  strontium.  The  half-power  points  §*   for  the  saturation  curve  (re- 


lated to  E  ,  see  Chapter  2)  is  inversely  proportional  to  the  quantum 


X 

r  qpp  Than for  7\     k  inuofcol 

efficiency  of  a  flame  or  plasma.  For  a  given  element  and  excitation 
line,  here  strontium  at  460.7  nm,  the  ICP  was  expected  to  have  a  higher 
quantum  efficiency  than  the  flame.  This  was  due  to  the  prevailing  argon 
atmosphere  of  the  ICP,  as  opposed  to  the  predominantly  molecular  atmo- 
sphere of  the  flame.  The  saturation  curves  (Figure  11)  indicated  that 
the  quantum  efficiency  of  strontium  in  the  ICP  was  slightly  higher  than  in  the 
flame.  A  second  measure  of  the  quantum  efficiency  of  strontium  in  the 
flame  or  plasma  was  the  saturation  broadened  excitation  profile  (Appen- 
dix 1  and  Reference  (62)).  A  broader  profile  indicates  a  higher  quantum 
efficiency,  as  shown  by: 

p(U  1/2 

SA    =  6\.  (1  + — ) 

exc    L  v    s,.  J 
p  UJ 


where 


5A    =  full  width  at  half  maximum  (FWHM)  of  the  excitation  profile 


exc 


of  the  atom  (nm) 


6X.  =  FWHM  of  the  Lorentzian  (undisturbed)  profile  of  the  atom 
(nm) 
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o 
p(?0  =  power  density  of  the  laser  (W/m  ) 

s  ? 

p  (X  )  =  saturation  power  density  (W/m  ) 

The  term  p  (x  )  is  inversely  proportional  to  the  quantum  efficiency  (see 

Chapter  2).  Therefore,  the  6X   depends  upon  quantum  efficiency. 

Saturation  broadened  SX   are  given  below  for  strontium  in  the  flame 

and  in  the  ICP  under  identical  power  density  conditions. 

Atom  Cell  ^exc  (full  power) 

flame  0.66 

ICP  1.13 

The  increase  in  <SX    in  the  ICP  also  bears  out  the  expected  increase  in 
quantum  efficiency  in  the  argon  atmosphere. 

Hydrogen  Based  Flames 

Hydrogen  based,  rather  than  hydrocarbon  based,  flames  were  chosen 
for  the  major  investigation  of  saturation  because  of  the  expected  im- 
provement in  saturation  due  to  a  much  less  complex  and  reactive  small 
molecule  (quenching)  mixture  in  the  flame  gas.  Two  inert  gas  diluents 
were  investigated  to  find  the  role  they  play  in  quantum  efficiency  in 
flames.  Nitrogen  was  chosen  for  its  obvious  place  in  any  air  oxidant 
flame.  Nitrogen  is  a  good  diatomic  quencher.  It  is  often  replaced  for 
atomic  fluorescence  purposes  by  a  monatomic  species,  argon,  which  is  a 
very  poor  excited  state  quencher.  Therefore,  it  allows  a  higher  quantum 
efficiency.  This  yields  a  higher  sensitivity  in  application  of  atomic 
fluorescence  to  analysis. 

Saturation  curves  for  the  two  diluent  gases  are  shown  in  Figure  12. 
The  two  curves  are  very  close  in  plateau  values  and  in  shape.  This  was 
an  indication  that  the  change  in  flame  gas  had  very  little  effect  in 
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either  the  atomization  of  strontium  or  the  quenching  of  excited  state 
strontium  atoms.  Relative  information  was  derived  for  Y  and  ny  from 
this  plot.  Relative  saturation  powers  (proportional  to  E^)  were  found 

A 

from  the  points  at  which  the  fluorescence  was  50%  of  the  maximum.  The 
saturation  powers  are  proportional  to  1/Y  (see  Chapter  2).  In  addition, 
the  slopes  of  the  1/Bp  vs.  1/E^  curves  were  taken  as  measures  of  1/Y  also. 
The  values  for  R„s  and  Ry£  were  computed,  and  are  shown  below. 
RYS  1.2 

RYE  1.1 

The  plateaus  of  fluorescence  were  used  as  relative  measures  of  nT 
(see  Chapter  2).  The  intercept  values  of  the  1/Bp  vs.  1/E  curves  were 
used  also  as  relative  measures  of  n-p  The  values  for  R...  and  RNB  are 
shown  below. 

RNI  °'91 

RNB  0.87 

The  agreement  of  the  two  methods  for  finding  relative  Y's  and  n's 

for  the  two  flames  was  quite  good.  The  percent  relative  error  for  the 

Y  measurement  was  9%.  The  percent  relative  error  for  n_  was  less  than 

J/o  • 

The  errors  in  the  ratios  are  probably  random  in  nature.  The  percent 
relative  standard  deviation  for  these  measurements  was  about  5-10%. 

Absolute  Y's  were  measured  by  means  of  the  E?  method  using  equation 
12  from  Chapter  2  under  the  conditions  of  Figure  12,  and  additionally 
under  the  conditions  of  Figure  13  using  the  flame-shielded  flame  and 
the  small  slit  mentioned  in  Appendix  1  placed  in  front  of  the  photodiode. 
The  results  are  shown  below.  Both  flame  stoichiometrics  were  the 
same. 
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time 


Figure  14.  Composite  boxcar  scan  of  fluorescence  pulse  of 

strontium.  The  gate  was  halted  at  several  delay 

times  (DT)  to  record  the  saturation  curves  seen 

in  Figure  15.  X  =  X       =   460.7  nm. 
ex   em 
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Strontium  Y~,  ,in  capillary  burner         0.50 
Strontium  Y?,  ,in  flame-shielded  flame      0.24 
As  the  results  indicate,  careful  attention  must  be  placed  on  select- 
ing a  correct  and  representative  aperture.  The  capillary  burner,  with 
no  flame  shield,  certainly  had  an  edge  effect,  producing  a  substantially 
poorer  quantum  efficiency.  This  was  entirely  masked,  however,  by  the 
inadequacy  of  the  measurement  of  the  E^  caused  by  using  an  aperture 
larger  than  the  laser  beam.  As  can  be  seen  in  equation  12,  the  quantum 
efficiency  is  inversely  proportional  to  E,  ,  and  therefore  directly  pro- 
portional to  the  measured  laser  beam  cross-sectional  area.  The  approxi- 
mately 2-fold  increase  of  quantum  efficiency  corresponded  in  this  case 
to  the  same  decrease  in  the  estimate  of  E:. 

A 

Total  number  density  was  figured  from  the  slope  of  the  1/Bp  vs. 

1/E  curve  after  calibration  of  the  intercept,  and  then  the  slope,  under 

12 
the  conditions  of  Figure  13.  This  procedure  resulted  in  n_  =  1.3  x  10 

_3 
cm  ,  a  very  reasonable  nT  for  100  ppm  Sr  in  a  hydrogen  based  flame.  A 

second  value,  derived  from  equation  10,  yielded  a  value  of  n,  =  3.1  x  10 

_3 
cm  .  This  value  agreed  very  well  with  the  prior  value.  The  small 

discrepancy  could  have  been  due  to  calibration  of  the  boxcar  plug-in 

module  or  the  photometer.  The  boxcar  module  had  been  under  repair  prior 

to  these  experiments,  but  had  not  been  calibrated.  Unfortunately  the 

equipment  required  for  a  check  of  this  value  was  not  available. 

A  temporal  scan  of  the  fluorescence  pulse  from  the  photomultipl ier 

produced  an  interesting  observation  for  strontium.  Several  overlapping 

scans  were  taken  to  define  the  composite  temporal  shape  shown  in 

Figure  14.  The  delay  time  for  the  boxcar  was  set  at  4  different  points 

on  the  peak,  and  4  sets  of  saturation  measurements  were  taken.  These 
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data  are  shown  in  Figure  15.  As  expected,  the  leading  edge  of  the  pulse 
was  not  saturated.  Surprisingly,  however,  the  trailing  edge  of  the  pulse 
was  saturated.  Normally,  assuming  two  levels  available  for  laser  and 
thermal  processes,  a  decrease  in  signal  could  not  be  observed  while  the 
system  was  saturated.  This  meant  that  a  third,  presumably  photon- 
induced,  process  was  necessary  to  explain  the  phenomenon.  Two  possi- 
bilities existed  for  the  explanation  of  these  data,  compound  formation 
and  laser-induced  ionization.  Compound  formation  was  rejected  as  a  cause 
of  this  effect,  as  a  two  body  collision  would  be  necessary.  Presumably 
the  only  suitable  species  available  in  this  flame  for  compound  formation 
would  be  oxygen  atoms  or  hydroxyl  molecules.  Even  with  favorable  orien- 
tation of  the  collisions,  they  probably  could  not  happen  fast  enough  to 
account  for  the  subnanosecond  phenomena  which  were  seen  here.  A  more 
probable  explanation  was  the  absorption  of  a  second  photon  of  wavelength 
460.7  nm,  which  could  bring  the  atom  within  0.3  eV  of  its  ionization 
potential.  A  thermal  event,  not  so  selective  as  a  compound  forming 
collision,  could  then  easily  ionize  the  atom.  Unfortunately,  an  experi- 
ment to  measure  electron  current  generated  via  laser-induced  ionization 
with  the  aid  of  the  nitrogen  pumped  dye  laser  was  unsuccessful  (63). 
In  private  communication  with  the  investigator  it  was  found  that  poor 
results  were  attributed  to  laser  pulse  width  considerations.  (While 
the  instantaneous  power  produced  ions,  the  time  spread  on  the  electron 
bunch  produced  only  a  small  signal.  The  best  type  of  laser  for  this 
experiment  has  been  found  to  be  one  with  a  longer  pulse  width  and 
highpower,  i.e.,  high  pulse  energy.) 
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S odium 

Measurement  of  sodium  parameters  was  hampered  by  an  experimental 

inconsistency.  In  all  other  measurements,  the  monochromator  bandwidth 

was  narrow  enough  to  exclude  any  other  radiating  wavelengths  from  the 

saturated  atoms.  For  sodium,  however,  both  resonance  lines  fell  within 

the  monochromator  bandpass.  Estimation  of  the  n,  values  from  BF   and 

max 
continuum  absorption  required  inclusion  of  both  sodium  resonance  lines 

in  the  theoretical  treatment.  The  nT  value  for  sodium  was  arrived  at 

via  two  different  routes.  A  value  was  calculated  from  the  saturation 

plateau  for  the  argon  diluted  flame  by  the  use  of  equation  10  from 

Chapter  2.  In  addition,  a  continuum  absorption  experiment  was  used  to 

give  an  independent  measure  of  ru,  as  described  by  DeGalan  and  Winefordner 

(17).  The  lack  of  a  suitable  bandwidth  monochromator  (ax    fc  5  x  6\   .  ) 

v  mono       atom' 

in  a  suitable  position  necessitated  substitution  of  an  argon  ion  pumped 
dye  laser  as  a  continuum  source  for  this  technique.  The  dye  laser  out- 
put bandwidth  was  measured  by  the  use  of  a  high-dispersion  monochromator 
(used  in  the  second  order).  The  optical  bandwidth  was  0.174  nm.  As 
the  sodium  atomic  absorption  linewidth  was  only  ca.  0.005  nm  (61),  this 
laser  was  considered  a  continuum  source.  The  monochromator  bandwidth 
described  in  (17)  was  taken  as  the  laser  bandwidth.  Measurements  of  a 
(fraction  absorbed)  were  made  in  both  the  argon  and  nitrogen  diluted 
flames.  In  addition,  saturation  curves  were  taken  in  both  flames 
(Figure  16).  The  n  's  calculated  from  both  methods  are  shown  below. 

Method  Value  (1  ppm) 

10    ? 
continuum  absorption  3.6  x  10   cm 

Bp  2.9  x  1011  cm"3 
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The  agreement  between  these  values  was  considered  quite  good  in 
light  of  the  fact  that  the  methods  are  independent,  and  several  calibra- 
tions of  different  types  are  involved.  Calibrations  of  absolute  optical 
radiances  are  considered  to  be  quite  difficult  and  error  prone.  An 
agreement  of  2-fold  is  considered  state  of  the  art  for  inexperienced 
laboratories. 

The  value  for  R...  was  determined,  and  augmented  by  an  RN  derived 
from  the  continuum  measurements.  However,  in  the  case  of  sodium  in  the 
nitrogen  diluted  flame,  the  saturation  was  incomplete.  Therefore,  the 
R^B  value  could  not  be  calculated.  Ratio  values  (R«)  are  given  below. 
The  values  were  corrected  for  differing  aspiration  rates  for  the  two 
f  1  ames . 

Method  Value 

RN,  continuum  absorption  0.89 

R„,  intercepts  0.84 

The  agreement  is  quite  good. 

Since  the  ratios  calculated  by  these  two  methods  did  agree  so  well, 
the  intercept  ratio  was  used  to  calculate  the  saturation  plateau  value 
for  the  nitrogen  diluted  flame.  This  enabled  evaluation  of  the  E? 
parameter,  which  was  otherwise  unavailable.  This  allowed  a  second 
means  of  determining  RY.  The  values  are  shown  below. 

Rye  5.3 

RYE  14.0 

Calcium 

In  Figure  17,  saturation  curves  for  calcium  in  HOA  and  HON  flames 
are  given.  The  Ry  values  for  the  two  flames  via  two  methods  are  given 
below. 
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Rys  1.31 

RYE  1.83 

A  Y  value  was  calculated  from  equation  12  for  the  HOA  flame,  and  is 
shown  below. 

YH0AE  °-52 

An  n,.  value  calculated  from  equation  10  is  shown  below. 

nTH0AB  2'6  x  1qU  cm"3 

Values  of  R..  derived  from  the  two  methods  are  shown  below. 

RNI  0.81 

RNB  °-83 


Indium 

Indium  was  marginally  saturated  in  the  stoichiometric  air-acetylene 
flame  (Figure  18),  which  is  a  highly  quenching  flame  with  respect  to  the 
HOA  flame. 

In  contrast  with  this  flame  Figure  19  shows  the  saturation  curve  in 
the  HOA  flame.  The  much  lower  power  and  much  more  complete  saturation 
shown  in  the  resonance  ( A  ex  =  xem  =  410  nm)  measurement  in  this  figure 
indicated  a  much  higher  quantum  efficiency.  In  addition,  the  observation 
of  x32  (451  nm)  fluorescence  followed  the  expected  behavior  in  showing 
the  same  saturation  shape  as  A...  observation.  The  2  ■*■  3  transition 
(451  nm)  could  also  be  marginally  saturated,  as  seen  in  the  same  figure. 
The  measurement  of  power  for  these  curves,  however,  was  not  reliable 
enough  to  provide  anything  but  the  relative  data  given  above. 

A  more  reliable  power  measurement  was  made  for  the  data  shown  in 
Figure  20.  The  values  for  Y  and  n-j.  in  the  HOA  flame  given  below  were 
arrived  at  by  the  calibration  of  the  1/BF  vs.  1/E  curve. 
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YHOA  °'28 

nTH0A  2.0  x  1010  cm"3  (1  yg/ml ) 

An  Rv  value  was  available  only  from  the  slope  method.  The  plateaus  of 
the  two  flames  were  not  developed  enough  to  find  Es  values. 

A 

Rys  1.03 

An  R^  value  was  only  available  from  the  intercept  method,  and  is  shown 
below. 

RNI  1.05 


Noise  Power  Density 

1/2 
Noise  power  density  (W/Hz  '  )  measurements  preliminary  to  satura- 
tion noise  power  density  were  made  under  the  conditions  shown  in  Figures 
21-26.  Since  the  measurement  of  saturation  noise  power  density  will 
only  be  possible  in  the  lower  frequency  ranges  (up  to  about  50  Hz,  due 
to  repetition  rate  limitations  on  pulsed  lasers),  the  low  frequency 
range  (0-50  Hz)  data  collected  were  chosen  for  these  figures.  The 
higher  collection  rate  (higher  frequency  spread)  data  had  no  interesting 
features  except  for  harmonics  of  60  Hz,  and  a  band  of  discrete  fre- 
quencies (centered  on  6  kHz)  introduced  through  the  electronics  power 
supply. 
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CHAPTER  5 
RESULTS  AND  DISCUSSION 

Quantum  Efficiency  and  Total  Number  Density 

Quantum  efficiency,  Y,  and  total  number  density,  nT,  values  are 
notably  lacking  from  the  literature,  though  they  are  of  great  interest 
to  analytical  chemists,  physicists,  and  engineers.  The  reason  for  this 
lack  of  literature  support  is  the  tedium  of  this  type  of  measurement, 
as  mentioned  in  Chapter  1.  Nevertheless,  several  items  are  available 
for  comparisons  in  flames  of  the  type  investigated  here. 

The  quantum  efficiency  (or  yield  factor)  has  been  previously  esti- 
mated in  flames  yery   similar  to  the  HOA  and  HON  flames  used  here  (64). A  value 
for  Y.  for  Sr  in  an  HOA  flame  yery   close  in  proportions  to  the  one  in 
this  work  was  0.16,  as  compared  to  0.24  in  this  work.  In  addition,  the 
ratio  of  this  Y  to  a  Y  evaluated  in  a  HON  flame  (gas  proportions  some- 
what different  from  this  work)  was  2.3,  as  compared  with  1.2  from  this 
work.  The  agreement  here  is  very   good,  considering  that  the  initial 
use  of  laser  excited  saturated  atomic  fluorescence  for  these  measurements 
was  in  this  work.  The  HOA  flame  Y  values  are  remarkably  similar.  The 
difference  in  the  ratio  values  can  be  explained  by  the  increased  pro- 
portion of  the  strong  quencher  nitrogen  in  the  flame  used  in  Reference 
(64). 

Another  source  of  literature  data  of  this  type  is  the  work  by 
Hooymayers  and  Alkemade  (65).  Data  on  sodium  quantum  efficiency  in  HO 
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type  flames  were  measured  by  these  authors.  When  values  for  Y  were  taken 
for  the  two  flames  with  gas  ratios  closest  to  those  in  this  work,  the 
RY  value  was  8.2.  The  Ry  values  found  via  two  different  methods  in 
this  work  were  Rys  =  5.3  and  RyE  =  14.  The  numerical  average  of  these 
is  9.5,  which  is  quite  close  to  the  published  value.  The  ratio  derived 
from  E^  for  the  two  flames,  14,  is  based  on  an  indirect  calculation  of 

nTH0N  from  nTH0A'  as  the  Plateau  condition  was  not  reached  in  the  HON 
flame.  Therefore,  the  RyE  is  somewhat  suspect.  The  Rys,  5.3,  is  some- 
what more  substantial,  and  does  not  differ  significantly  from  the  value 
in  Reference  (65).  As  the  quantum  efficiency  is  a  strong  function  of 
gas  composition  (even  with  the  same  diluent  gas),  the  difference  between 
the  slope  data  ratio  and  the  published  data  was  quite  reasonable. 

The  absolute  value  of  the  published  YHQA  (65)  for  Na is  0.32.  The  value 
arrived  at  from  E^  in  the  present  work  is  0.03.  If  the  estimation  of 
power  were  in  error  for  any  reason,  this  would  be  reflected  in  the  Y 
value  found  by  this  method. 

The  quantum  efficiency  for  calcium  in  the  HOA  flame  from  the  pre- 
sent work,  0.52,  seemed  quite  reasonable,  in  light  of  the  quantum 
efficiency  seen  for  strontium,  which  should  behave  similarly  in  the 
flame.  No  literature  data  were  available  for  direct  comparison, 
however. 

The  quantum  efficiency  for  indium  was  of  the  same  order  of  magni- 
tude as  strontium  and  calcium,  and  therefore  seemed  reasonable.  Again, 
no  direct  literature  comparison  could  be  made. 

The  Ry  values  seemed  to  be  rather  lower,  in  general,  than  expected. 
Nitrogen  was  expected  to  be  a  much  better  quencher  than  argon.  In  addi- 
tion, the  concentrations  of  the  two  diluent  gases  were  predominant  in  the 
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flarnes,  leading  to  the  presumption  that  they  would  dominate  the  quenching 
rate.  In  the  sodium  case,  however,  the  Rv  value  was  several  times 
higher,  more  in  line  with  expectations.  The  literature  available  per- 
tains to  sodium  (64)  and  strontium  (65).  Quenching  cross-sections 
(a,  A  )  in  HO  type  flames  are  given  below  for  both  atoms  in  Table  5. 


Table  5 

Quenching  Cross-Sections  for  Sodium  and  Strontium 
in  Hydrogen  Flames 


Na  (64)  Sr  (65) 

Quencher  a  (A)2  a   (A)2 


N2  40  21 

02  62  400 

H20                           1  67 

Ar                           2.3  <  1 


The  sodium  case  (high  R  )  seemed  to  fit  the  expectation  that  a 
large  quenching  change  could  be  expected  from  a  change  of  the  major 
diluent  (<?Ar/aN  =   17).  The  other  main  constituent  of  the  burnt  gases, 
water,  had  a  low  a  (1  A  ),  and,  therefore,  did  not  affect  the  quenching 
significantly. 

In  the  strontium  case  (R  =   1);  the  a   ratio  is  still  high.  Water, 
however,  has  a  larger  cross-section  (ca.  3-fold  higher)  for  strontium 
than  nitrogen.  Therefore,  the  influence  of  the  change  of  diluent  gas 
will  be  much  less  than  in  the  case  of  sodium.  In  addition,  the  cross- 
section  of  oxygen  was  20-fold  higher  than  that  of  nitrogen.  Since  ^ery 
small  changes  in  the  oxidant  to  fuel  ratio  around  the  stoichiometric 
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mixture  used  here  can  grossly  affect  oxygen  concentration,  this  could 
greatly  affect  the  trend  caused  by  nitrogen. 

The  RN  values  for  this  work  were  all  quite  close  to  1,  indicating  a 
very  small  role  for  nitrogen  in  the  atomization  process.  It  would  be 
expected  that  a  large  change  in  the  flame  gas  composition  from  a 
monatomic  constituent  (argon)  to  a  diatomic  diluent  (nitrogen)  should 
make  a  gross  difference  in  the  atomization  process,  through  thermal 
equilibration.  Apparently  this  does  not  happen  in  these  flames,  how- 
ever. Here  again,  no  comparative  data  were  available  in  the  literature 
for  the  two  diluent  gases. 

Noise  Power  Density 

A  very  preliminary  evaluation  of  the  results  of  the  noise  power 
experiments  is  given  below  in  Table  6.  Experimental  conditions  for  the 
data  are  found  in  Figures  21-26.  Changes  in  different  types  of  noise  are 
indicated  by  arrows  pointing  in  the  direction  of  change. 

Table  6 

Noise  Changes  Induced  by  Adding  an  Inert  Gas  Sheath 


Wavelength  Region 


Noise  Types 


1/f         White       Proportional 


OH  (306  nm)  J  |  \ 

C2  (516.5  nm)  —  t  t 

Sr  (460.7  nm)  —  f  f 
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The  OH  region,  as  was  expected,  exhibited  a  lower  white  noise  com- 
ponent when  an  inert  gas  sheath  was  added.  The  inert  sheath  minimized 
entrapment  of   02  at  the  edge  of  the  flame,  producing  less  OH.  This,  in 
turn,  caused  the  white  noise  (shot  noise)  component  to  decrease.  In- 
creased turbulence  caused  the  1/f  noise  and  also  the  proportional  noise 
to  increase.  Proportional  (or  whistle)  noise  was  nonexistent  in  any 
case  until  the  sheath  was  added.  It  appeared  in  all  cases  when  the 
sheath  was  added. 

The  C2  region  was  enhanced  by  addition  of  the  sheath,  since  a  richer 
atmosphere  favored  C,,  production.  The  white  noise  therefore  increased. 
The  low  frequency  (1/f)  noise  remained  the  same,  since  turbulence  at 
the  flame  edge  did  not  grossly  affect  C?  production. 

Strontium  1/f  noise  remained  the  same  with  addition  of  the  sheath. 
The  lack  of  a  change  in  1/f  noise  could  have  been  due  to  a  cancellation 
of  two  effects.  The  overall  strontium  signal  was  3-fold  lower  when  the 
sheath  was  on.  This  could  have  resulted  in  lower  1/f  noise,  which  would 
have  been  cancelled  by  an  increase  due  to  increased  turbulence. 


CHAPTER  6 
CONCLUSIONS  AND  FURTHER  WORK 

The  use  of  a  nitrogen -pumped  dye  laser  was  investigated  for  the 
measurement  of  flame  diagnostic  parameters.  The  assurance  that  all  the 
conditions  necessary  for  these  measurements  were  fulfilled  was  a  very 
time-consuming  and  tedious  project.  The  following  points  were  especially 
critical  in  development  of  the  method. 

System  Temporal  Response 

The  establishment  of  atomic  fluorescence  pulse  widths  via  reliable 
means  (streak  camera)  enabled  us  to  choose  the  correct  measurement 
system  (75  ps  aperture)  from  the  choices  available.  The  existence  of  a 
steady  state  temporal  plateau,  while  certain  to  occur  somewhere  within 
the  pulse,  had  to  be  long  enough  to  measure  without  averaging  unsaturated 
portions  of  the  pulse.  The  photomultiplier  rise  time  had  to  be  fast 
enough  to  follow  the  pulse  without  distorting  its  shape  unduly. 

Complete  Saturation 

Complete  saturation  was  desirable  for  several  reasons.  It  eliminated 

concern  over  the  necessity  for  spatial  homogeneity  in  illumination  of  the 

flame  volume  under  investigation.  Complete  saturation  also  provided 

BF   and  a  way  to  find  E?  for  easy  calculation  of  nT  and  Y.  These  values 
rmax  A 


•105- 


■106- 


(Bp   and  E,)  also  provided  confirmatory  data  for  ratios  of  n_  and  Y 

max     '  ' 

in  the  different  flames. 


Measurement  of  Ex 

Laser  power,  bandwidth,  and  area  also  had  to  be  evaluated  very 
carefully  in  order  to  provide  realistic  absolute  data  for  Y.  The  latter 
two  quantities  do  not  affect  ratios  of  Y,  of  course.  As  seen  in  the 
strontium  section  of  Chapter  5,  an  improper  evaluation  of  E,  can  easily 

A 

cause  a  3-fold  change  in  absolute  Y  estimation.  Laser  spectral  band- 

idth  had  to  be  evaluated  in  order  to  satisfy  the  requirement  for  a 
spectral  continuum  for  illumination,  also. 

Photometer  Calibration 


w 


The  photometer  had  to  be  calibrated  in  order  to  find  Bp    in  an 

max 


absolute  way,  to  get  nT. 

General  Considerations 

Careful  corrections  had  to  be  made  for  laser  scatter  into  the 
monochromator  at  high  irradiances.  Typically, scatter  from  the  laser 
contributed  significantly  to  the  observed  signal  until  laser  irradiance 
was  reduced  by  about  2  orders  of  magnitude. 

The  data  evaluated,  after  these  parameters  were  controlled,  were 
reasonable  and  generally  in  agreement  with  literature  values  of  nT  and 
Y  available. 

A  complication  arose  when  the  discovery  of  a  further  photoinduced 
process  came  about.  Photoionization  and  thermally  assisted  photoinduced 
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ionization  were  discussed  as  possible  explanations  for  these  data.  Two 
photons  of  the  same  wavelength  (460.7  ran)  were  enough  to  bring  strontium 
close  to  ionization  (0.3  eV  from  the  ionization  potential).  For  indium 
two  photons  (410  nm)  were  enough  to  exceed  the  ionization  potential. 
Further  work  by  another  student  found  recombination  ion  lines  upon 
scanning  the  emission  monochromator  while  illuminating  indium  with  the 
laser.  The  current  theory  does  not  include  the  processes  of  ionization 
or  recombination.  At  laser  powers  of  the  order  necessary  for  saturation 
of  some  elements  these  effects  may  have  to  be  considered  in  the  future. 

The  possibility  exists,  with  the  advent  of  the  excimer  laser,  that 
peak  powers  can  be  had  at  visible  and  ultraviolet  wavelengths  to  com- 
pletely ionize  atoms  in  a  reasonable  volume  in  a  flame.  If  this  power, 
and  the  longer  pulse  widths  of  this  type  of  laser,  can  completely  and 
selectively  ionize  a  given  atomic  population,  good  sensitivity  and  a 
certain  immunity  to  kinetic  (recombination)  factors  should  result.  This 
type  of  saturation  should  readily  avail  itself  of  theoretical  calibra- 
tion via  geometric  collection  efficiency  factors. 

A  second  possibility  for  measurement  of  nT  in  three-level  atoms  is 
to  use  two  dye  lasers  to  saturate  both  transitions  and  to  thereby 
equalize  populations  (according  to  their  degeneracies).  This  would 
allow  measurement  with  no  assumptions  about  kinetic  quenching  constants. 

Noise  Power  Density 

Since  only  preliminary  noise  power  spectra  have  been  measured  so 
far,  much  work  remains.  The  expected  reduction  in  noise  during  satura- 
tion has  yet  to  be  observed.  Further  characterization  remains  to  be 
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done  on  the  temporal  behavior  of  the  different  atoms  under  laser  excita- 
tion. It  is  hoped  that  the  application  of  saturation  to  noise  power 
spectra  measurement  will  allow  real  progress  to  come  about  in  the  critical 
areas  of  nebulization  and  burner  design. 

General  Comments 

Saturated  atomic  fluorescence  has  been  developed  and  critically 
evaluated  as  a  tool  for  flame  diagnostic  measurements.  It  promises  to 
be  more  straightforward  and  less  tedious  than  other  methods  of  arriving 
at  values  for  Y  and  rip  Hopefully  saturated  atomic  fluorescence  will 
enable  the  evaluation  of  noises  basic  to  design  features  in  atomic 
analytical  instrumentation. 

Spatial  and  temporal  information  obtainable  through  the  means  of 
saturated  atomic  fluorescence  are  available  by  no  other  means. 


APPENDIX  1 
SUBSIDIARY  MEASUREMENTS 

Since  the  measurement  system  used  in  these  studies  was  to  be  used 
for  the  measurement  of  fast  signals,  it  was  necessary  to  characterize 
its  response  on  a  very  short  (ns)  time  scale.  The  temporal  response 
of  the  system  was  investigated  in  several  ways,  using  several  auxiliary 
systems  to  the  one  mentioned  previously.  The  points  which  needed  to 
be  investigated  were 

1.  PMT  pulse  response, 

2.  laser  pulse  width, 

3.  atomic  fluorescence  pulse  width. 

Photomultiplier  Pulse  Response 

The  response  of  the  photomultiplier  tube  and  base  combination 
was  measured  by  the  use  of  a  small  laboratory-built  nitrogen-pumped 
dye  laser.  The  laser  pulse  duration  had  previously  been  measured  at 
about  600  ps.  This  pulse  was  directed  through  the  same  collection 
system  used  for  the  saturation  measurements.  For  all  practical  pur- 
poses, this  represented  an  impulse  of  light  to  the  photomultiplier 
tube.  The  pulse  time  response  of  the  PMT  was  then  measured  by 
scanning  a  75  ps  gate  across  the  output  pulse.  The  result,  shown 
in  Figure  27,  was  a  pulse  response  width  of  2.9  ns. 
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Figure  27.  Boxcar  temporal  scan  of  photomultiplier  pulse 

response.  Laser  input  pulse  width  was  ca.  600  ps, 
Photomultiplier  response  width  was  2.9  ns. 
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Laser  Pulse  Width 

The  temporal  features  of  the  laser  output  were  checked  in  two 
independent  ways.  The  output  of  a  calibrated  vacuum  photodiode 
(risetime  I  500  ps)  was  scanned  using  the  75  ps  gate  (Figure 28).  The 
data  collected  in  this  manner  were  checked  against  the  manufacturers 
calibration  via  the  use  of  a  borrowed  joulemeter.  The  joulemeter 
response  indicating  pulse  energy  could  be  reconciled  with  peak  power 
measurements,  made  via  connection  of  the  output  of  the  photodiode 
through  a  50  Q   load  to  a  fast  oscilloscope  (Type  454,  Tektronix, 
Portland,  OR).  The  pulse  energy  and  peak  power  could  be  related  if  a 
pulse  shape  and  pulse  width  were  found.  Pulse  width  from  the  photo- 
diode-boxcar  measurements  did  not  agree  with  the  manufacturer's 
calibration  data,  using  the  peak  energy  found  with  the  joulemeter. 
However,  using  a  variety  of  other  sources,  pulsed  and  continuous 
wave,  the  diode  calibration  was  upheld.  Therefore,  a  need  to  further 
investigate  the  laser  temporal  distribution  was  expressed,  as  the  laser 
pulse  width  and  time  response  of  the  PMT  to  the  fluorescence  caused  by 
the  pulse  needed  to  be  understood  thoroughly. 

Manufacturer's  literature  (66)  and  several  publications  featuring 
the  use  of  this  laser  list  the  duration  of  the  dye  laser  pulse  width 
as  4-5  ns.  The  joulemeter  data  indicated  times  ranging  from  300  ps 
to  approximately  2  ns  as  pulse  widths  for  different  dyes.  As  these 
were  considered  somewhat  nonsensical  for  this  laser  system,  it  was 
decided  to  use  a  very  fast  response  system  for  time  dependence  of  the 
dye  laser  pulses.  A  Hamamatsu  streak  camera  (Temporal  Dispersor, 
Hamamatsu,  Middlesex,  NJ)  was  used  in  conjunction  with  a  microprocessor 
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Figure  28.  Boxcar  temporal  scan  of  photodiode  pulse  response. 
Laser  input  pulse  width  was  ca.  600  ps.  Photodiode 
response  width  was  2.2  ns . 
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based  data  collection  and  display  system  (Temporal  Analyzer,  Hamamatsu,  Mid- 
dlesex ,-  NJ) .  This  system  could  trigger  the  laser  and  the  streak  camera 
with  the  use  of  suitable  delay  lines  and  trigger  generators  to  ensure 
overlap  of  the  camera  firing  and  the  laser  pulse.  The  laser  was 
directed  through  suitable  neutral  density  filters  (to  avoid  damage  to 
the  very  sensitive  streak  tube)  and  into  a  monochromator  via  two 
front  surface  mirrors.  The  exit  slit  of  the  monochromator  was  focused 
by  streak  camera  optics  on   the  streak  camera  entrance  slit.  Several 
streaks  were  observed  to  find  the  general  temporal  shape  for  the  given 
dye.  As  these  experiments  were  done  during  the  daytime,  severe 
triggering  jitter  existed.  Only  10%  or  so  of  the  streaks  collected 
were  usable.  Several  representative  streaks  were  collected.  Each  of 
these  was  recorded  via  a  digital  to  analog  converter  and  a  strip  chart 
recorder  (Model  SRG,  Sargent-Welch,  Cleveland,   OH).  The  results 
for  several  dyes  are  shown  in  Figures  29-32.  The  results  indicated 
that  the  dye  laser  system  used  here  was  either  producing  much  shorter 
and  more  temporally  complicated  pulses  than  expected,  or  a  serious 
systematic  error  was  being  committed.  Therefore,  the  temporal  output 
of  a  second  dye  laser  (DL2000,  Lambda-Physik,  Interactive  Radiation, 
Northvale,  NJ),  an  oscillator-amplifier  configuration,  was  also  in- 
vestigated. This  laser  was  expected  to  give  a  temporal  pulse  width 
in  the  same  4-5  ns  range.  The  streak  results  are  shown  in  Figure  33. 
This  result  was  as  expected,  and  the  other  results  were  therefore 
trusted.  The  calibration  of  the  photodiode  via  the  joulemeter  could 
then  be  reconciled  with  the  manufacturer's  calibration.* 


*Note:  Further  data  collected  since  experimental  work  for  this  disser- 
tation was  finished  indicate  that  variations  in  pulse  width  and  jitter 
of  the  laser  depend  strongly  on  repetition  rate  and  gas  flow  velocity. 
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Figure  29.  Streak  of  Molectron  dye  laser  scatter  (7D4MC) , 
1.5  ns  per  inch. 
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Figure  30.     Streak  of  Molectrondye  laser  scatter  (DPS),  1.5  ns 
per  inch. 
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Figure  31.  Streak  of  Molectron  dye  laser  scatter  (R6G)  ,  1.5  ns 
per  inch. 
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Figure  32.  Streak  of  Molectron  dye  laser  scatter  (bis-MSB) 
1.5  ns  per  inch. 
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Completeness  of  Saturation 

The  completeness  of  saturation  was  shown  for  the  several  atomic 
species  chosen  by  the  plateau  region  in  the  plot  of  fluorescence 
radiance  vs.  laser  irradiance  (Figures  9-20).  In  addition,  this 
provided  proof  that  the  precautions  against  spatial  laser  beam  in- 
homogeneity  were  working. 

Laser  Power  Measurement 

Laser  power  (peak  power)  was  measured  by  the  use  of  the  above- 
mentioned  vacuum  photodiode.  This  was  placed  in  front  of  the  laser 
with  a  large  enough  calibrated  neutral  density  filter  to  provide  a 
signal  still  linear  with  respect  to  the  power  (typically  4.0  neutral 
density).  The  output  was  connected  through  a  50  n  attenuator,  and  the 
peak  voltage  was  read  on  the  oscilloscope.  The  laser  power  could  be 
arrived  at  after  consultation  of  the  manufacturer's  calibration  table. 

Laser  Bandwidth 

The  laser  bandwidth  calibration  was  approached  in  two  separate 
ways.  When  the  first  bandwidth  data  were  taken,  a  suitable  mono- 
chrcmator  (high  resolution)  was  not  accessible  for  direct  measurement 
of  the  bandwidth  of  the  laser  output.  Inasmuch  as  saturation  of  an 
atomic  transition  can  give  a  larger  fluorescence  excitation  bandwidth 
than  the  naturally  occurring  bandwidth,  an  experiment  was  undertaken 
to  find  the  excitation  bandwidth  of  the  laser.  After  a  transition 
was  found,  a  wavelength  scan  of  the  laser  output  was  made.  The 
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Figure  33.     Streak  of  Lambda- Phys i k  dye  laser  scatter   (7D4MC), 
8.7  ns  per  inch. 
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fluorescence,  collected  by  the  same  low-resolution  monochromator  was 
that  of  the  integrated  emission  bandwidth.  A  suitable  theory  was 
developed  in  Omenetto  et  al .  (62)  to  understand  this  manifestation  of 
saturation  broadening  in  terms  of  the  laser  wavelength  spread  and  laser 
power  density.  The  low  power  limit  of  the  excitation  band- 
width vs.  laser  power  phenomena  turns  out  to  be  the  laser  bandwidth, 
as  might  be  surmised.  Plots  of  excitation  bandwidth  low  and  high 
power  limits  are  shown  for  the  dyes  used  in  this  study  (Figures  34-37). 

In  order  to  confirm  the  low  power  limits  as  a  source  of  laser 
bandwidths,  a  suitable  monochromator  (HR1000,  Jobin-Yvon,  Instruments 
SA,  Inc.,  Metuchen,  NJ)  was  used.  Resolution  (bandwidth)  of  the  mono- 
chromator was  found  by  wavelength  scanning  a  mercury  line  from  a  low 
pressure  pen  lamp  source  in  the  second  order.  This  resulted  in  a 

o 

0.12  A  bandwidth  for  the  monochromator.  The  dye  laser  outputs  for 
the  pertinent  atomic  transitions  were  scanned  directly  (with  suitable 
neutral  density  filters).  These  scans  were  done  with  the  use  of  the 
scanning  attachment  for  the  Molectron  dye  laser,  as  this  was  easier 
than  monochromator  scanning.  The  bandwidth  results  for  both  measure- 
ments are  shown  in  Table  7.*  In  several  years  of  use  the  anti- 
reflective  coating  of  the  dye  laser  cells  had  apparently  been 
destroyed.  This  allowed  the  development  of  a  second  laser  cavity, 
coupled  to  the  front  surface  of  the  dye  cell,  which  resonated  at  a 
different  wavelength,  and  with  a  different  temporal  position.  This 


*Note:  An  anomalous  side-lobe  structure  is  seen  in  the  results  of  both 
monochromator  scan  and  saturation  broadening  bandwidth  for  the  dye 
7D4MC.  This  structure  is  probably  related  to  the  anomalous  second 
pulse  seen  in  some  temporal  scans  of  this  dye. 
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Figure  34.     Excitation  broadening  of  strontium.     Relative  atomic 

fluorescence   (vertical   axis)   vs.  wavelength   (horizontal 
axis),  obtained  by  scanning  the  dye  laser  wavelength  past 
the  atomic  wavelength   (460.7   nm) ,  while  continuously 
recording  the  atomic  fluorescence  signal    (7D4MC  dye). 

a.  low  power  limit—represents   laser  bandwidth. 

b.  high  power  limit—represents   broadened 
sensitivity  under  saturation  conditions. 
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Figure  35.     Excitation  broadening  of  indium.     Relative  atomic 

fluorescence   (vertical   axis)   vs.  wavelength   (horizontal 
axis)    (DPS  dye,  center  wavelength,  410  nm). 

a.  low  power  limit—laser  bandwidth. 

b.  high  power  limit. 
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Figure    36.     Excitation  broadening  of  sodium.     Relative  atomic 

fluorescence   (vertical   axis)   vs.  wavelength   (horizontal 
axis)    (R6G  dye,  center  wavelength,  589.0  run). 

a.  high  power  limit. 

b.  low  power  limit—laser  bandwidth. 
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Figure  37.  Excitation  broadening  of  calcium.  Relative  atomic 

fluorescence  (vertical  axis)  vs.  wavelength  (horizontal 
axis)  (bis-MSB  dye,  center  wavelength,  422.7  nm) . 

a.  low  power  limit— laser  bandwidth 

b.  high  power  limit 
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Table  7 
Excitation  Profile  vs.  Scanned  Laser  Bandwidth 


sx  (A) 


Element         Direct         Fluorescence  Excitation  Profile 
Measurement      — 


Ar/02/H2  N2/02/H, 

Ca  0.23  0.20  0.20 

Sr  0.23  0.42  0.23 

Na  0.36  0.23  0.40 

In  0.24  0.19  0.28 
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observation  was  further  born  out  by  the  abbreviated  pulse  width  of  the 
other  dyes,  and  the  limiting  spectral  bandwidths  exhibited.  These  were 
approximately  three  times  those  specified  by  the  manufacturer.  The 
further  evidence  of  some  temporal  mode  structures  on  most  of  the 
dyes'  output  seemed  to  concur  with  the  above.  The  times  between  maxima 
of  the  temporal  modes  corresponded  to  a  distance  of  about  1  cm,  the 
width  of  the  dye  cell.  This  indicated  a  strong  coupling  within  the 
cell  with  an  additional  coupling  to  the  entire  cavity.  Further  ex- 
periments will  be  done  with  a  new  dye  cell  to  confirm  this  line  of 
reasoning.  In  the  meantime,  the  sensitivity  of  the  temporal  width  to 
tuning  of  the  cavity  invites  the  goal  of  very  short  laser  pulses  from 
a  commercial  system  not  designed  for  this  type  of  work.  In  atomic 
spectroscopy,  this  invites  the  study  of  rate  constants  of  quenching, 
etc. 

Laser  Cross-Section 

The  evaluation  of  the  cross -sect ion  involved  in  the  excitation 
of  atoms  in  the  viewing  area  was  purely  based  on  the  geometry  of 
collection.  The  vertical  distance  was  the  slit  width,  as  the  mono- 
chromator  was  fixed  on  its  side.  The  horizontal  distance  was  assumed 
to  be  the  width  of  the  inner  (seeded)  flame,  which  was  somewhat  smaller 
than  the  laser  beam.  Therefore  a  mask  was  placed  over  the  face  of  the 
photodiode  before  peak  power  was  measured.  The  aperture  measured 
0.5  mm  by  4  mm,  and  was  centered  in  the  laser  beam,  which  had  been 
centered  at  the  correct  vertical  and  horizontal  displacement  during 
alignment.  The  homogeneity  of  the  laser  beam  is  discussed  in  the 
following  section. 
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Laser  Seam  Homogeneity 

The  necessity  of  a  spatially  homogeneous  laser  beam  for  meaningful 
diagnostic  measurements  has  been  discussed  at  length  previously  (67). 
In  order  to  use  the  data  from  these  experiments,  it  was  necessary  to 
prove  the  homogeneity  via  direct  sampling,  as  well  as  by  the  results 
of  the  diagnostic  experiments.  The  laser  beam,  suitably  filtered  as 
previously  described  was  measured  at  the  edge  of  the  flame  cell  farthest 
from  the  laser.  The  same  vacuum  photodiode  as  used  previously  was 
masked  with  a  100  ym  pinhole  and  suitable  neutral  density  filters. 
It  was  then  translated  through  the  beam  in  two  dimensions  via  a 
laboratory-built  screw  translator  mechanism.  The  output  of  the 
photodiode  was  measured  with  the  boxcar  system  using  the  164  plug-in. 
The  resulting  voltage  was  recorded  on  a  strip  chart  recorder.  These 
data  were  then  entered  manually  into  a  digital  computer.  The  data 
were  plotted  using  a  three-dimensional  plotting  routine  developed  by 
the  author  (see  Appendix  2).  These  results  are  shown  in  Figure  38. 
The  homogeneity  in  the  central  portion  was  quite  good.  This  central 
portion  (0.5  mm  by  4  mm)  was  the  portion  directly  exciting  the  field 
of  view. 

Self-Absorption 

Self-absorption  was  proved  to  be  negligible  for  all  measure- 
ments by  introducing  a  10-fold  higher  concentration  of  analyte  into 
the  flame  with  a  1.0  neutral  density  filter  inserted  between  the 
flame  and  the  PMT.  In  all  cases,  this  produced  the  same  size  signal 
as  the  10-fold  lower  concentration.  The  lower  concentration  was  used, 
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however,  for  the  measurements  because  of  the  tendency  for  higher  con- 
centrations to  reduce  the  degree  of  saturation  (68). 

Spectral  Continuum 

The  laser  bandwidth  of  approximately  0.25  &   was  assumed  to  be  a 
continuum  despite  the  possibility  of  a  saturation  broadened  absorption 
profile.  Power  broadening  should  not  be  present  at  these  laser  power 
levels. 

Achievement  of  Steady  State 

The  achievement  of  a  steady  state  in  the  population  of  the  ex- 
cited state  is  necessary  to  measure  a  meaningful  quantum  efficiency 
or  total  number  density.  While  steady  state  could  be  achieved 
in  at  least  two  ways,  only  one  of  these  was  suitable  for  these 
measurements.  The  first  steady  state  has  been  predicted  (67)  as  a 
condition  in  which  (in  a  three-level  system)  the  populations  of  the 
ground  state  and  the  highest  excited  state  are  equalized  by  the  fast 
rising  edge  of  a  laser  pulse.  This  steady  state  relaxes  to  a  second 
steady  state  plateau  involving  the  second  level.  It  is  this  feature 
which  accurately  reflects  the  quantities  that  we  are  interested 
in.  As  a  means  of  determining  that  such  a  steady  state  exists,  a 
temporal  scan  (streak)  was  made  of  the  fluorescence  of  the  highest 
excited  state  of  indium  as  a  function  of  laser  excitation  tuning. 
The  tuning,  as  mentioned  before,  strongly  affected  the  pulse  width 
of  the  dye  used.  Fluorescence  corresponding  to  X.-  excitation 
and  a23  excitation  are  shown  in  Figures  39  and  40.  In  order  to 
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Figure  39.  Streak  of  indium  atomic  fluorescence.  Excitation 

wavelength  410  nm  (see  streak  in  Figure  41).  Emission 
wavelength  451  nm.  Concentration  1000  ug/ml.  19  ns 
per  inch. 
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Figure  40.  Streak  of  indium  atomic  fluorescence.  Excitation 

wavelength  451  nm  (see  streak  in  Figure  33).  Emission 
wavelength  410  nm.  Concentration  1000  yg/ml.  19  ns 
per  inch. 
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Figure  41.  Streak  of  Molectron  dye  laser  scatter  (DPS),  8.7  ns  per 
inch. 
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collect  these  streaks  a  very   high  (ca.  1000  ug/ml)  concentration 
solution  was  necessary.  In  addition,  the  non-resonance  transition 
had  to  be  used,  as  scatter  was  too  high  for  resonance  measurements. 
This  transition  also  had  to  be  in  steady  state  for  the  measurements 
to  be  meaningful.  No  evidence  of  peaking  in  the  highest  excited  state 
was  found.  It  was  evident  that  the  temporal  width  of  the  fluorescence 
was  far  longer  than  that  of  the  laser,  and  that  a  significant  plateau 
existed  at  the  top  of  the  pulse.  The  steady  state  was  quite  evident. 
It  was  also  evident  that  either  the  threshhold  of  saturation  had  to 
be  quite  low  on  the  laser  excitation  pulse  or  that  the  combined 
decay  rate  connecting  the  third  excited  state  with  the  other  states 
had  to  be  very   low.  If  the  latter  were  true,  the  plateau  would  be 
much  shorter  than  it  is,  as  the  laser  drops  off  with  quite  a  fast 
fall  time  compared  to  the  length  of  the  plateau  (see  Figure  41). 
Therefore,  it  can  be  considered  that  saturation  occurred  quite  early 
in  the  pulse.  Although  there  is  no  data  representing  the  relative 
timing  of  the  leading  edges  of  the  laser  and  the  fluorescence  pulses, 
it  could  be  assumed  that  if  the  laser  shape  were  not  changed  grossly 
that  the  falling  edge  of  the  atomic  fluorescence  pulse  represented  a 
relaxation  of  the  excited  state  population.  A  later  section  on  rate 
constants  discusses  this  possibility. 

Coherence  Effects 

Coherence  effects  (population  changes  synchronous  with  the 
electric  field  of  the  laser  output)  would  destroy  the  approach  used 
in  the  theoretical  section  to  find  quantum  efficiency  and  total  number 
density.  These  effects  are  not  expected  to  be  seen  in  high  pressure 
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environments  at  the  power  levels  used  here,  as  collisions  of  the 
dephasing  type  are  quite  frequent  (47). 

Measurement  of  Fluorescence  Depth 

Since  the  laser  beam  exceeded  the  width  of  the  inner  (seeded) 
flame,  the  fluorescence  collection  depth  was  taken  as  the  width  of  the 
flame  in  the  collection  axis,  4  mm. 

.Calibration  of  Photometer 

A  DC  calibration  procedure  was  used  for  calibration  of  the 
photometer  used  for  collection  of  pulsed  data.  Pulsed  calibration 
procedures  and  calibrated  sources  are  not  available.  A  tungsten  strip 
lamp  (GE  pyrometer  6818,  S6,  GE,  Nela  Park  Cleveland,  OH)  calibrated 
against  a  National  Bureau  of  Standards  calibrated  strip  lamp  (EPUV 
1068,  National  Bureau  of  Standards,  Washington,  D.C.)  was  placed  at 
the  focus  of  the  optical  system  for  fluorescence  collection.  The 
lamp  was  operated  at  the  standard  35  A  current  (AC)  specified,  and 
allowed  to  stabilize  for  1/2  hour,  shielded  on  all  sides  by  drop 
cloths.  Measurement  of  output  current  was  made  under  the  experimental 
conditions  at  the  wavelengths  specified  in  the  lamp  calibration  table. 
The  resulting  data  are  shown  in  Figure  42. 

Quenching  Rate  Constants 

The  atomic  fluorescence  streak  data  exhibit  the  plateau  mentioned 
in  the  Steady  State  section  of  this  appendix.  This  feature  means  that 
if  the  saturating  pulse  is  short  enough,  yet  intense  enough  to  populate 
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substantially  the  highest  state,  that  a  combined  quenching  and 
radiative  lifetime  for  that  level  could  be  found  from  streak  data. 
The  signal  to  noise  ratio  of  the  present  data  is  quite  poor,  but  if 
a  smoothed  decay  profile  were  used,  it  could  establish  an  upper  limit 
for  a  decay  lifetime.  The  problem  with  this  approach  is  that  the 
laser  temporal  response  cannot  be  seen  with  enough  amplitude  resolu- 
tion (with  the  present  equipment)  to  determine  if  some  response  to 
the  tail  of  the  pulse  is  determining  the  observed  decay.  In  addi- 
tion, with  the  jitter  problems  as  they  stand,  no  temporal  relationship 
could  be  established  between  onset  of  lasing  and  onset  of  fluor- 
escence, although  they  should  coincide  closely.  A  computer  acqui- 
sition of  several  temporal  streaks  of  the  laser,  and  then  the  fluor- 
escence, followed  by  alignment  of  the  leading  edges  has  been  tried 
by  one  author  (69)  to  eliminate  problems  of  this  sort.  In  the 
present  case,  this  was  out  of  the  question,  as  the  jitter  was  far 
too  large. 

Noise  Power  Density  Measurement 

The  system  used  to  collect  these  data  has  been  used  before  (69), 
but  will  be  described  here  briefly.  Figure  43  depicts  the  system. 
The  source  of  the  noise  for  the  spectra  collected  was  an  ordinary 
laboratory  flame  (nebulizer  and  spray  chamber:  Perkin-Elmer  Corpora- 
tion, Norwalk,  Conn.)  with  a  laboratory-built  capillary  burner 
(Haraguchi  et  al .  (70))  and  an  inert  gas  sheath  (laboratory-built). 
A  stabilized  gas  flow  system  similar  to  the  one  mentioned  previously 
was  used.  Emission  signals  were  collected  by  a  simple  point  to  point 
lens  system  at  the  slit  of  a  medium  resolution  monochromator 


Figure  43.  Noise  power  density  collection  system.  Component  symbols 
are  as  follows:  F  =  flame,  L  =  lens,  MO  =  monochromator , 
PMT  =  photomultiplier  tube,  PS  =  power  supply,  R  =  load 
resistor,  ACA  =  AC  amplifier,  C  =  blocking  capacitor, 
A/D  =  analog  to  digital  converter,  D/A  =  digital  to 
analog  converter,  MC  =  minicomputer,  0  =  oscilloscope, 
X-Y  =  X-Y  recorder,  T  =  terminal. 
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(EU-700,  McPherson  Instrument  Co.,  Acton,  Ma.).  A  1P28  PMT  (RCA, 
ElectroOptics  and  Devices  Div.,  Lancaster,  Pa.)  was  used  at  -700  V 
(Model  240,  High  Voltage  Supply,  Keithley,  Cleveland,  OH).  The 
output  current  was  converted  to  a  voltage  with  one  of  two  load 
resistors  (1  kft  or  16.2  kn).  An  AC  amplifier  was  used  to  block  the 
DC  voltage  and  amplify  the  AC  or  noise  signal,  in  order  to  best  use 
the  dynamic  range  of  the  data  acquisition  system.  The  output  of  this 
amplifier  was  taken,  after  suitable  adjustment  of  an  upper  cut-off 
RC  filter,  through  another  blocking  capacitor  (1  pF,  mylar,  used  to 
block  DC  offset  from  the  amplifier)  to  an  analog  to  digital  converter 
(LPS-11,  Digital  Equipment  Corporation,  Maynard,  Mass.).  This  ADC 
was  interrogated  by  the  computer  at  programmable  intervals  to  collect 
2048  samples  of  noise  voltage.  These  were  Fourier  transformed.  The 
average  of  100  Fourier  transforms  at  a  given  acquisition  rate  were 
converted  to  a  power  spectrum,  which  could  be  displayed  in  several 
ways  (see  Figure  43).  The  digital  data  were  transferred  in  a 
parallel  fashion  to  another  digital  computer  (PDP-11/20,  Digital 
Equipment  Corporation).  These  data  were  stored  in  binary  fashion 
for  fast  acquisition  by  a  sequential  data  plotting  program  (see 
Appendix  2). 


APPENDIX  2 
COMPUTER  PROGRAMMING  AND  PLOTTING 

A  useful  and  necessary  function  of  a  computer  system  is  the  ability 
to  express  the  results  of  its  computations  in  a  way  which  is  intel- 
ligible to  the  programmer  and  the  user.  It  has  been  said  that  man, 
since  he  lives  in  a  three-dimensional  world,  has  a  limited  ability  to 
visualize  in  higher  dimensions.  Mechanical  considerations  limit  his 
ability  to  display  data  in  more  than  two.  At  best,  he  can  produce 
a  pseudo-three-dimensional  display. 

The  PDP  11/20  system  used  for  some  of  this  research  has  lacked 
for  many  years  even  the  most  rudimentary  display  system.  A  need  was 
established  for  plotted  data  from  several  sources,  of  two  types, 
sequential  (or  channel)  type  data,  and  paired  X-Y  data.  Both  types  lend 
themselves  to  two  and  (pseudo-)three-dimensional  plotting.  Therefore, 
a  project  was  undertaken  to  produce  a  stepping  motor  plotter  to  ful- 
fill the  need  for  this  type  of  display  and  for  creation  of  simple 
diagrams. 

An  old  X-Y  analog  plotter  chassis  was  scavenged  for  the  bed  and 
pen  movement.  Stepping  motors  were  attached  via  gears  and  belts. 
Driving  cards  were  installed,  as  well  as  end  of  travel  micro- 
switches. 

As  our  computer  had  no  general  digital  input-output  function,  a 
request  was  made  to  local  software  support  specalists  to  install  and/or 
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modify  existing  routines  as  part  of  the  group's  BASIC  language  package. 
This  included  a  stepper  drive  command  capable  of  addressing  eight 
separate  motors.  This  enabled  the  plotter  to  be  programmed  in  a 
somewhat  wasteful  way  from  BASIC.  A  modular  program  of  providing 
plotting  packages  was  commenced  at  that  time.  The  basis  of  the  total 
package  was  a  45  line  subroutine  which  plotted  a  best  line  between  two 
points  on  the  plotter.  No  provision  was  made  for  runoff,  so  scaling 
and  manipulation  of  data  had  to  provide  for  the  plotter  size  internal 
to  a  driving  program  utilizing  the  best  line  routine.  The  routine 
could  have  been  made  shorter,  but  included  the  following  points  neces- 
sitated by  the  present  form  of  BASIC: 

1.  software  integerization  to  provide  symmetry  in  the  function, 

2.  skip-out  to  delta  X-delta  Y  change  when  the  pen  is  up, 

3.  skip-out  to  single  axis  call  when  only  one  axis  is  changed, 

4.  mini -vector  approximation  with  correction  after  each 
mini-vector  in  only  one  dimension. 

The  basic  best-line  routine  has  been  used  as  the  core  of  several 
data-plotting  routines: 

1.  sequential— e.g.,  optical  multichannel  analyzer,  noise 
power  spectra,  or  streak  camera  data, 

2.  X-Y  format--e.g.,  scans  of  intensity  vs.  a  single  dimension, 
as  in  profiling  cuts  of  a  flame  or  a  plasma, 

3.  X-Y-Z  format— e.g.,  scans  of  intensity  vs.  two  dimensions, 
as  in  multidimensional  profiling,  wavelength-intensity-time 
profiles,  etc. 

These  programs  were  used  on  various  types  of  computer  files,  as  suited 
each  program,  and  type  and  volume  of  data. 
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In  addition,  the  best-line  routine  served  the  etch-a-sketch  or 
box  diagram  routine  via  a  very  simple  language  containing  a  number  of 
symbols  for  commonly  used  figures  in  spectroscopic  diagrams.  This 
was  useful  for  slide  and  figure  preparation  for  talks  and  journal 
articles. 

A  third  utilization  of  this  routine  was  for  the  character  genera- 
tion module.  This  type  of  programming  has  generally  been  done  in 
assembly  language  to  facilitate  speed.  In  this  case,  the  generation 
has  all  been  done  in  BASIC  with  approximately  a  ten  times  reduction 
in  speed  from  a  commercial  system.  All  of  the  letters  and  numbers 
and  several  of  the  ASCII  symbols  have  been  provided  for  labelling  of 
charts,  etc.,  created  by  the  other  programs.  In  some  cases,  the 
character  generation  module  has  been  integrated  with  other  graphics 
modules. 

Specifics  follow  for  the  best-line  routine  and  programs  utilizing 

it. 

Best-Line  Routine 

Lines  5000-5020 

This  is  a  set-up  section.  It  saves  current  position  (C  is  the  X 
axis  position,  D  is  the  Y  axis  position)  in  CI  and  Dl.  Flags  J  and  K 
are  set  to  indicate  the  sign  of  change  for  each  axis,  for  the  purpose 
of  the  software  integerization. 

Line  5025 

This  line  transfers  to  a  short  form  when  the  pen  is  raised  (i.e., 
when  B(0)  =  1).  The  short  form  (lines  5900-5915)  steps  off  the  X  axis 
change  first,  then  the  Y  axis  change,  to  save  time. 
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Lines  5028,  5029 

These  lines  transfer  control  to  short  form  sections  when  only  one 
axis  is  changed,  again  to  save  time. 

Lines  5520-5530 

The  Y  axis  only  is  changed. 

Lines  5500-5520 

The  X  axis  only  is  changed. 

Lines  5030-5050 

These  lines  decide  whether  the  greater  change  is  in  the  X  axis 
or  the  Y  axis  and  steers  the  program  to  the  correct  subroutine.  In 
addition,  a  true  slope  is  computed  and  the  best  integer  approximation 
(absolute  value)  is  found. 

Lines  5600-5665 

This  is  the  best  line  generator  for  a  change  in  X  axis  larger 
than  a  change  in  Y  axis. 

Lines  5800-5870 

This  subroutine  generates  the  best  line  for  the  opposite  case, 
change  in  the  Y  axis  larger  than  the  change  in  the  X  axis.  The  in- 
verse of  the  original  slope  is  computed,  and  the  integer  of  the 
absolute  value  is  found.  The  integerized  slope  is  a  minivector,  when 
used  in  conjunction  with  the  direction  flags  J  and  K.  This  minivector 
is  the  closest  approximation  to  the  real  line  that  can  be  made  by 
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discrete  steps.  At  every  repetition,  it  is  correct  in  the  dimension 
on  the  bottom  of  the  slope  (in  this  case  the  X  axis  change)  since 
steps  of  unit  size  are  always  taken  on  this  axis.  On  the  other  axis, 
if  the  correct  position  is  kept  track  of,  a  correction  can  be  made 
whenever  the  correction  amounts  to  at  least  one  step  (as  computed  in 
lines  5830-5850).  If  the  minivector  is  repeated  as  many  times  as  there 
are  steps  between  the  points  on  the  lower  axis  (X  in  this  case),  the 
best  line  has  been  plotted.  This  is  the  basis  of  all  the  graphics 
functions.  New  position  is  provided  to  the  subroutine  in  XI  and  Yl. 
The  change  in  terms  of  steps  is  computed  from  the  present  position, 
C  and  D.  A  variable  list  can  be  found  in  Table  8. 

Sequential  Data  Plots 

0P1.BAS  Lines  1-4  and  1000-2000 

A  virtual  data  file  must  be  provided  with  values  ranging  from  0 
to  35,565.  Scaling  is  done  at  the  keyboard  to  provide  for  expansion. 

The  file  is  searched  for  the  maximum  integer  value. 

Line  1100 

The  actual  working  value  is  assigned  to  this  integer. 

Line  1105 

The  expansion  is  decided  here  by  the  actual  largest  value  to  be 
displayed.  All  values  topping  this  value  are  limited  to  it  (cropped) 
in  a  subroutine  (lines  1600-1610).  The  X  axis  is  treated  similarly, 
except  no  cropping  is  provided.  Therefore,  care  must  be  taken  when 
expanding  via  lines  1070-1090. 
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Table  8 
Variable  List  for  Best-Line  Routine 


CI  Saves  present  X  position 

Dl  Saves  present  Y  position 

C  Present  X  position 

D  Present  Y  position 

J  Delta  X  sign  flag 

K  Delta  Y  sign  flag 

XI  New  X  position 

Yl  New  Y  position 

B  Digital  output  matrix  for  word  1 

51  Actual  slope 

52  Integerized  absolute  slope 
P9  Loop  variable  for  best-line 

56  Actual  pen  position 

57  Difference  between  actual  and  correct  pen  positions 
M  Sign  flag  for  S7  for  X/Y 

L  Sign  flag  for  S7  for  Y/X 
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Lines  1110-1390 

These  lines  provide  for  axis  drawing  and  tic  marks  on  either  axis 
or  both.  The  position  of  the  tic  marks  is  remembered  for  later 
labelling. 

Lines  1400-1500 

These  lines  provide  for  actual  plotting,  using  the  best-line 
routine  and  either  placing  dots  or  drawing  lines  between  each  data 
point. 

Axis  Labelling  and  Numbering 

0P1.BAS  Lines  200-760 

Lines  200-351 

These  lines  allow  numbering  at  remembered  tic  mark  positions  and 
multiple  line  labellingfor  the  X  axis. 

Lines  355-465 

The  same  function  is  performed  for  the  Y  axis. 

Lines  500-760 

This  is  the  caption  section  subroutine,  which  services  the 
labeling  and  numbering  sections.  It  cuts  up  the  caption  into  in- 
dividual ASCII  characters,  finds  out  the  size  of  the  caption,  and 
steps  back  half  of  this  length  to  center  it.  The  captions  can  be 
plotted  horizontally  or  vertically,  and  character  size  is  at  the 
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operator's  choice.     This  routine  makes  use  of  the  best-line  routine 
via  the  character  generation  routine. 

Character  Generation  Subroutine 

Lines  4400-4456 

This  is  the  look-up  table  which  provides  the  look-up  in  variable 
J8  for  the  character  mini  vectors.     This  variable  accesses   a  binary 
virtual  data  file  containing  the  mini  vectors   and  pen  up-down  commands 
comprising  the  57  ASCII   characters   and  special   symbols.     These  symbols 
were  coded  and  converted  into  the  binary  virtual   file  via  the  program 
LEDAT.BAS.     The  file  name  is   LEI. DAT.     The  data  file  may  be  output  in 
an  intelligible  manner  by  running  VFP1.BAS,  and  referring  to  the 
character  sort  section  of  this  subroutine. 

The  manner  in  which  the  virtual   data  file  is  used  is  straight- 
forward.    The  file  is  arranged  as   follows: 

Item  1— the  number  of  symbols  contained  in  the  file. 
Item  2— the  number  of  symbols  plus  one--the  beginning  number 
of  the  minivector  data  for  the  corresponding  symbol. 
Items  through  the  end  of  the  file-symbol  minivector  tables  as 
follows: 

the  first  item  in  each  table  is  the  number  of  vectors 
in  the  symbol,  starting  at  the  lower  left  corner  of  a 
7  by  7  dot  matrix,  and  ending  at  the  lower  right 
corner.     The  following  items  are  minivector  com- 
ponents by  threes, 
1.     Pen  up   (0)   or  down   (1)    bit 


-148- 


2.  X  change  in  steps 

3.  Y  change  in  steps. 

The  character  routine  uses  the  BASIC  table  to  find  the  virtual 
file  element  pointing  to  the  minivector  set.  A  loop  is  run  through 
the  correct  number  of  times  (specified  by  the  first  element),  each 
time  using  the  pen  set  routine  and  scaling  the  minivector  as  desired 
before  calling  the  best-line  section  to  draft  the  pen  stroke.  At  the 
finish  of  the  letter,  the  pen  is  moved  a  space  corresponding  to  half 
a  letter.  Then,  the  loop  is  used  for  the  next  letter. 

Data  Passage  Between  Computers 

Passage  of  sequential  data  files  from  the  PDP  11/34  in  a  parallel 
fashion  to  the  PDP  11/20  was  accomplished  via  the  PSTRN  program  (11/34) 
and  the  REC  program  (11/20).  The  only  formality  necessary  for  com- 
munication in  this  fashion  is  the  handshake  routine.  This  was  per- 
formed on  the  11/34  end  with  "peek"  and  "poke"  commands.  The  11/20 
used  the  augmented  BASIC  commands  "INPB"  and  "PLSB."  A  1000  by  16 
bit  word  transfer  takes  about  1/2  minute,  which  is  not  ^/ery   spectac- 
ular, but  quite  easy  to  program.  The  major  holdup  was  in  BASIC, 
which  was  quite  slow.  An  assembly  language  program  transferring  to 
the  11/34  from  a  hardware  device  ran  about  an  order  of  magnitude 
faster.  The  ease  of  use  of  this  method  balanced  the  time  consumed. 
The  transfer  of  ten  or  more  files  at  a  stretch  would  be  accomplished 
with  no  intervention. 

Passage  the  other  way  has  not  been  developed  in  great  detail 
yet.  Test  programs  have  run  with  only  hardware-caused  errors. 
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Plans  are  under  way  for  bit-parallel   character-serial   passage 
of  ASCII   files  as  a  completely  general  method  of  exchange  for  programs 
and  data  files.     BASIC   (on  the  11/20)   has  ASCII   coding  and  decoding 
statements   for  alphanumeric  sending  and  receiving.      Fortran   (on  the 
11/34)   does  not  have  these  features.     However,  a  look-up  table 
should  suffice. 
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Best-line  routine--  This  subroutine  serves  all  plotting 

functions . 


■>  5  0  y  0  ) 


STEP" (K 


0) 


5000  C1=C\D1=IJ 

5005  J=l 

5010  IFr  (X.l-C)>-~OGO  10  5015  \J~-J 

5015  K=l 

5020  IF  (Y1-D)>=0G0  TO  5025  \K~-1 

5025  IF    B<0)=1G0  TO  5020 

5026  GOSUB  5900  \RETURN 

5028  IF    X1=CG0  TO  5520 

5029  IF  Y1=DG0  10  5500 

5030  S  J .  =  (  X 1-C  )  /  (  Y  1.  -D  )  \S2=INT  <  ADO  (  0  J.  )  • 
5040  I!-  ABS(Sl)>=iGO  TO  5600 
5050  GOSUB  5000  \G0  TO  5400 
5400  RETURN 

5500  CALL  "SIEP"  (INT  (ABS>(  XI -C)  !' ;■'■:  l»50»l) 
5505  C=C+IKT(ADS(X1-C) >*J 
5510  GO  TO  54  00 

5520  CALL  '  STEP "< INT  CABS ( Yl-B) ) 
5525  D=D+INT(ABS(Yl-rO  >*K 
5530  GO    TO  54  00 
5600  FOK    P9-1  TO  ADS(Yl-D) 
5610  CALL  ,STKP,(S2*J»S0»1)\CALI 
5625  C  =  C  +  S  2  *  J  \  D  -  D + K 
5630  S6=P9#ABS<S1  )*JiCl 
5640  S7=S6--C 

5644  M--1 

5645  IF  S7>--"0G0  TO  5650  \tf~-l 
5650  CALL  " SI LF  '  ( INT ( APS ( 07 )  > *M , 00  .  1  ) 
56^5  C=C+INT(ABS<S7) )*M 
5660  NEXT  P? 
5665  GO  TO  54  00  ' 
5300  S1  =  1/S1\S2"INT(AI:S(S1  )  ) 
5805  FOR  P9--1  TO  ABS(Xl-C) 
5810  CALL  'STEP" (S2&Kr 50*0) \C 
5820  C«C+J\D=D+S2#K 
5830  S6=P9*A!:<'S(SJ.  )Jf-.KiHi 

5840  S7=S6-B 

5044  L-l 

504  5  IF  S7>~000  TO  0050  \l...-  =  ~.l 

5050  CALL  ■STir-  (j'N'K  AI.:G<  07)  )*L?50,0'> 

5355  D~IH  it'j  (r,l'C(0/)  )  :■;  L 

5860  NEXT  P9 

5070  RETURN 

5900  CALL  ,SrEP'(IMT(rtCS(Xl-C))^Jr50»l)\C-CiINT<ADS(Xi-C)>*J 


')LL  "OTL!"'  (  J?  00 «  1  ) 
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0P1.BAS--A  section  of  the  program  dealing  with  sequential  plotting 


DEF 
DEF 

DEF 
GO 


1  DIM 

FNA(A'B»C)~A*BX:200/C 
FNB(W»X>  Y>Z)-=WX.X*Y/2 
FNC(D>ErF,GylO~D*E*F/(G 


Uc 


BETWEEN 
•iLUE  IN 
M.UE  ON 
SCALE  (IMC!-! 
Y  VALUE  IN 


X 

•!  IT 


EC) 


TO  J. OOO 

1000  PRINT  TILE  NAMED  FOR  INPU1 

1001  PRINT  * (UP  TO  1024  INTEGERS 
1005  A-0 

1010  PRINT  'FIEF  NAME'  J  \  INPUT  7?'!- 

1020  OPEN  79%    FOR  INPUT  AS  I  LI..  L 

1025  PRINT  '  t  OF  POT  Ml  C  TO  PLOT" 

1030  FOR  J6-0  TO  X3-1 

1040  ]  F  0  F  1  (  J  6  )  < = A  G  0  1  0  1  0  6  0 

1050  A=VF.1(J6) 

1060  NEXT  J 6 

1070  PRINT  '#  OF  STEP 

1080  PRINT  'MAXIMUM  X 

1090  PRINT  'MAXIMUM  X 

1095  PRINT  UY  AX I 

1100  PRINT  "MAXIM 

1105  PRINT 

1110  PRINT 

1120  IF  Y4$='N"G0 

1130  30 SUB  4  900 

1140  X1=FNB(X2» 

1145  Y1=0\G0SUB 

1150  Xl=-50 

1160  GO SUB  5000 

1170  Y1=FNA(Y2> Y6 7 Y6) 

CALL  rSTEP"(Ylf3C 

CALL  "STEP" (~Y1 -70? 

GOSUB  4950 

PRINT  "TIC  MARKS  (Y 

IF  T*--BN"GO  TO  14  00 

P=l 

PRINT  "X  VALUE" »\IN 


B  (  1  5  )  f  R  (  1  0  )  r  C  (  1  0  ) 

I) 

■it.' ST 


BE 

JES 

A    BA 

0    TO 

1.02 
JT 

V  ■  ! 

RTUAL  INTEGER  FJLI 

)  " 


MAXIMUM 
DRAW  AX!: 
TC 


:s.X6 

5000 


13  30 
1190 
1200 
1250 
1260 
1265 
1270 
1275 
1230 
12S5 
1290 
1  295 
1300 
1310 
1320 
1325 
1330 
1340 


Y  VAL 
S  (  Y 

14  00 

rX5) 

0  ) 


ON 

N  ) 


U: 


N) 


JT 


IF  X7--< 
X1=FNB(X2.X 
R(P)=X1 
GOSUB  5000 


:c 


TO 

X7 


X7 
.Y1----0 


GOSUB 

CALL 

GOSUB 

CALL  ' 

P=P+1 

GO  TO 

0=1 


4900 

STEP 
4950 
STEP' 


(•-20r30;0 


(20 


1270 


1350  PRINT 
1355  IF  Y7- 
1360  Y1--FN,' 
1365  S(Q)"Y1 
1370  GOSUB  50c 


VAU 

900 


IE" 
TO 


(Y2»Y 


0) 


:npu 
ioo 

.  i  -•■  -- 


r  y 


'01  NTS'  ;  VINPU 
CORRESPONDS 
S" »\INPUT  X6 
?\INPUT  Y2 
CORRESPONDS  TC 


Q"J\INPU1  X5 


.INPUT  Y 


vMNPUT 


INPUT  T 


o  \oosu; 


00 


-152- 


1375 
13  80 

1385 

Crl 

PL 
LL 
0  f  1 

"STEP 

"0  1  yj- 

* (-20, 
" (20f 3 

30  v  l 
Of  I) 

>\GG 

GUI 

;  19^ 

1390 

GO 

TO 

.1350 

1400 
1405 

1406 

XI 
PR 
IF 

INT 
L$  = 

ri=o\ 

"  L I N 

30SUB 
■3    OR 
J  TO  1. 

50  00 
DOTS 
4  1 0 

(1... 

")R 

D )  "  P 

1407 

14  or; 

1409 

Y.1 
CU 
GO 

-FNi 
3UB 

3un 

:(VF1  (0)  ,Y2 
1600  \oosu 

4  900 

•:i  50 

A.  •  Y6 
>'<0 

>  *v 

00 

1410 

r  u 

i  j<! 

■■-:'.     T( 

)  X3-1 

1420 

1430 

Yl 

-FNC 

-:J6* 

'(05  J 
:X2 

'  J6> , Y 

2f  Y5 

>■  fi »  Y, 

> )  >: 

2  0  0 

1435 
144  0 

GO 
GO! 

-U..li.: 

1600 

50 ''10 

1450 

IF 

L'l.- 

"L  "G( 

)  TO  i: 

500 

1460 
1500 

GOf 

ne; 

n  j 

4  9  0  0 
6 

\GOSM) 

' '.) 

1510 

00 i 

lUB 

4950 

1520 
1530 

J.  sj  O  vJ 

XI  = 

o  o  ;-. 

GLC 

=0\Y 
:UD 

OF 

j.  ■-  0 
5000 
Or  1 

1540 
1550 

PR  I 
IF 

i"  l" 
7  P  $ 

"RE  PL 

-  '  Y  '  G 

OT  (Y 
0  TO  1 

OR  !■ 

000 

)  " ;  \ 

IN 

"'UT  Z 

1560 
1570 
15S0 

IF 

GO 

NT 

7  0  :!; 

TO 

"CAPTION  ( 
='N'GO  TO  2 

10 

y  or 

000 

N )  " 

V  \ 

CNPUT 

1600 
1605 

1610 

IF 

y  i.  - 

■•>  L  T 

Y 1  < 
Y2* 
URN 

^Y2*2 
200 

DO  GO  T 

0  16 

.1.0 

2000 

END 

INPUT 
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0P1.BAS--A  subsection  of  this  program  dealing  with  axis 
labelling  and  generation  is  shown  here. 


X  AX  IB  AND  NUMBER 


■iXIS    CND    NUMBER 


200  P R INT  "  X  "  5 \ I N P U T  X 1 

210  PRINT  ,Y" ?\ INPUT  Yl 

220  GO SUB  5000 

230  RETURN 

250  PRINT  "CP ACE  BE THE CM 

255  Y1--U4 

260  FOR  Pl-1  TO  P-l 

270  X1=R(P1 ) 

230  GOSUB  5000 

290  GOSUB  500 

300  NEXT  P.1 

301  GO  TO  3:1.0 

305  PRINT  'ANOTHER  LINE  OF  X  AXIS  L, 

306  IF  V$O"Y"G0  TO  355 
310  PRINT  "SPACE  BETWEEN  X  AXIS  AND 
320  Yl=-W4 

330  Xl-FND(X??X3.X6?X5)/2 
340  GOSUB  5000 

350  GOSUB  500 

351  GO  TO  305 

355  PRINT  'SPACE  BETWEEN  Y 
360  Xl=-W4-50 
370  FOR  Pl-1  TO  Q-l 
380  Y1=S(P1) 
390  GOSUB  5000 
4  00  GOSUB  500 

410  NEXT  PI 

411  GO  TO  4  20 

415  PRINT  "ANSI  HER  LINE  OF 

416  IF  0$<>"Y"GO  TO  490 
420  PRINT  ' SPACE  BETWEEN 
430  Xl=-U4-50 
440  Yl=PNA(Y2?Y6vY6)/2 
450  GOSUB  5000 
460  GOSUB  500 
4  65  GO  TO  <U5 
490  Xl=0\Yi~0\ GOSUB  5000 
4  95  RETURN 

500  PRINT  "POSITION  P 
505  P6=0 
510  L-J  J  «C\W2~D 

520  PRINT  ■HORIZONTAL  (I!)  OR  {,' 
530  PRINT  \ PRINT  \ PRINT 

540  PRINT  "LET  I  EL  SIZE  (ABOUT  LOO  STEP 
550  PRINT  "CAPTION  *  ?\  INPUT  OK 


\ INPUT  W4 


.  (Y  OR  N) 
iXIS  LABEL 


J  YINf 

;\inf 


UT 

UT 


W4 


MNPU 


W4 


Y 


Y  AX 
AXIS 


;EL  (Y  OR  N) 


Y  AXIS  LABEL 


MNPL 
MNPU 


W4 


AT  CENTER  POINT  0 


(0 


LETTE 

:r  inc 


ON 


R I N 


MNP 


UT  LI 

cs 
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551  B4--0 

555  conuo  700 

560  IF  L.1*="V"UG  TO  5V0 

569  Y1=D 

570  Xl-C- ( 54/2 ) - ( C9/2 >  * ( LEN < CI  * ) -1 > /2 
575  GOSUB  5000 

580  GO  TO  600 

590  X1«C 

591  Y1~D~<  S4/2)  •-  (CB/2)  *  (LENCCl*  )  -1  '>  /2 
595  GOSUB  5000 

600  FOR  09=1  TO  I...  E  N  <  C  ].<■!; ) 

610  A2*--SEG*(C1*»0?»09> 

620  GOSUB  -4400 

650  NEXT  09 

640  X1«U1\Y1«W2\GGSUB  5000 

650  RETURN 

700  P6--1 

710  FOR  G9==l  TO  LEN<C1$) 

720  A2*--SEG$(C1?  r09>09) 

725  P6=l 

730  GOSUB  4400 

735  P6  =  C 

740  NEXT  09 

760  RETURN 
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0P1.BAS--A  subsection  of  the  program  dealing  with 
character  generation  is  shown  here. 


4400 

IF 

A2$- 

:>gA 

■GC 

TO 

4  4  0 1 

\JS 

=  1\G0 

";  0 

44  57 

4401 

IF 

A2*< 

:> "  p 

"GO 

TO 

4-102 

V.'O 

"2\G0 

TO 

445*7 

4402 

IF 

A2*< 

:>*c 

"  G  G 

TO 

4403 

\  J  0 

-3\G0 

i  0 

4  4  5  7 

4403 

IF 

A2*'; 

■•  '  D 

"GO 

TO 

•1  4  04 

\J8 

=4\G0 

TO 

4457 

4404 

IF 

A2*< 

>  "  E 

"GO 

TO 

4405 

\J0 

~5\G0 

TO 

4  4  5  /' 

4405 

IF 

A  2  $  ■•' 

;;,  "  F 

°  GO 

TO 

4  4  06 

\JG 

=6\G0 

TO 

44  57 

4406 

IF 

A2*< 

::■• "  G 

"GO 

TO 

A  4  ()  ~f 

\J8 

=  7\GC 

TO 

4457 

4407 

IF 

r-;  2  ■,>  *-'■ 

::•■■ "  hi 

'GO 

TO 

4400 

\JS 

-S\GO 

7  0 

4457 

4408 

IF 

A  2 '■!>■•■: 

>  "  I 

"GO 

TO 

440? 

\JS 

-9\G0 

TO 

44  57 

4409 

I  r 

A2J< 

••;.  »   1 

°  G  CJ 

TO 

4  410 

\J8 

=10\GC 

TO 

4457 

4410 

IF 

A  2  $  < 

>  "  K 

"GO 

TO 

4411 

\JS 

=  11\GG 

7  0 

4  4  5  7 

4411 

IF 

A  7  $  :- 

>  "  L 

•GO 

TO 

4412 

\J3 

=  12\G0 

ro 

4  4  57 

4  4 1  2 

IF 

A  2  -i>  •■'■ 

>  "  ri 

"GO 

TO 

4413 

\J3- 

"13\G0 

t  n 

44  57 

4413 

IF 

A2$< 

>  •  N 

"GO 

TO 

4414 

\J8= 

=14\G0 

TO 

44  57 

4414 

I F 

A2S< 

>  "  0 

'  GO 

TO 

4415 

\  JS- 

15XG0 

TO 

4457 

4415 

IF 

i-". '!'  *  •:'' 

> "  p 

'GO 

TO 

4  416 

\JS- 

■16\G0 

TO 

4457 

4416 

IF 

A  2$< 

;:.-■ "  Q 

GG 

TO 

4  4  1  7 

.\J3= 

=  17X00 

TO 

4457 

4417 

IF 

A2*< 

'•'• "  R 

GO 

TO 

4410 

\J8- 

=  13X00 

TO 

4  4  57 

44  IS 

IF 

A2i-< 

:-• "  S 

GO 

TO 

4  4  1  9 

\ , !  G  • 

:J.?\GO 

TO 

4457 

441? 

IF 

A  *"'  -' :- 

:■-  ■  T 

GO 

TO 

44  20 

\J3-- 

•20\GG 

TO 

4  457 

4420 

IF 

i'  i  2  '■!-  •'•. 

:  ■- "  U  ' 

GG 

TO 

4  421 

\JS= 

=  21\G0 

"0 

4  457 

4421 

I F 

A2iK 

> "  V ' 

GO 

TO 

4  422 

\JS- 

•22\G0 

)  0 

4457 

4422 

I F 

A  2  $  < 

:■- "  W  ' 

GO 

TO 

4  423 

SJS= 

23\G0 

TO 

4  4  57 

4423 

IF 

0  2  v  ■■'■'. 

:•- "  X " 

GO 

TO 

4  4  2  4 

\J8- 

24XG0 

ro 

44  57 

4424 

1 1" 

A2">< 

,  „  y  , 

GO 

TO 

4425 

\J&-- 

25\G0 

"1  0 

4  '■  :  / 

4425 

IF 

A  2  ^  <. 

:•  »  2  ' 

GO 

TO 

4  4  26 

\JS= 

26X00 

"0 

■1  157 

4426 

t  cr 

,*,  "-  -|.  ■■;' 

••■  "  1  " 

GO 

TO 

4  4  2  7 

\J3= 

27X00 

TO 

4457 

4427 

I F 

A  2  v  •■'" 

•  "  2 " 

GO 

TO 

4  4  20 

\Jfs- 

','    \i.  "i  i 

ro 

4  457 

4428 

I F 

A  2  !■•<: 

, .  3  . 

GO 

TO 

'  <  4  2  V 

\JS= 

2  V  \  1  1 0 

TO 

4457 

4  42? 

IF 

f  i  2  V1  r-. . 

:•  "  /,  " 

GO 

TO 

4130 

\JS~ 

30XGO 

TO 

4  4  57 

4  430 

]  r 

A2*<: 

•'•  "  Li  ' 

GO 

TO 

4  4  3 1 

\,JC  - 

3i\G0 

TO 

-1 4  5  7 

4431 

Ih 

A2$<' 

,   «  /,  -■ 

CO 

TO 

,<l  A   -•  "P 

\js- 

32XG0 

TO 

44  57 

-156- 


4432 

IF 

A2$< 

>"7"G0 

TO 

4433  .\JS=33\G0 

TO 

.4457 

4433 

IF 

A2$< 

>"8*G0 

TO 

4434  \J8=34\G0 

TO 

4  457 

4434 

J  F 

A2$< 

>"9"G0 

TO 

4435  \J3=35\G0 

TO 

44  57 

4435 

IF 

A2$< 

>"0"G0 

TO 

4436  \JG--=36\G0 

10 

4  4  57 

4436 

IF 

A2$< 

>" ( "GO 

TO 

4  137  \J8"37\G0 

TO 

4457 

4437 

IF 

A2$< 

:■■ " ) "GO 

TO 

4  4  30  \J8~-38\G0 

TO 

44  57 

443S 

IF 

A2$< 

::■"-«  GO 

TO 

443?  \J8-3V\G0 

TO 

4  457 

443? 

IF 

fa?*    :, 

:■"  .  "GO 

TO 

4440  \J3~40\G0 

TO 

4457 

4440 

IF 

A2$< 

:  "  »  "GO 

TO 

4441  \J8-41\G0 

TO 

44  57 

4441 

IF 

A2$< 

>'+"60 

TO 

4  44  2  \JG=42\G0 

TO 

445"; 

444  2 

IF 

A2$< 

>•/■ GO 

TO 

4  44  3  \J8=43\G0 

TO 

44G  ::' 

4443 

IF 

A2$-: 

>'3"G0 

TO 

4144  \J3=44\G0 

TO 

445  7 

44  44 

IF 

A2$< 

:"cl:  -go 

TO 

4445  \J8-=45\GC 

TO 

4  457 

4445 

IF 

A2$< 

::•"--•  "GO 

TO 

4-*46  \J8-46\G0 

ro 

4457 

4446 

If 

A2$< 

>">"G0 

TO 

44  47  \J8~47\G0 

TO 

4457 

4  447 

IF 

A2$< 

:•  "<"G0 

TO 

4  448  \JQ~43\G0 

TO 

4457 

444S 

I  r 

A2$< 

>"i'iE'GO 

TO 

4  4  4  9  \  J  8 ■-■  4 9  \  l'J  G 

TO 

44  57 

444? 

1 F 

A2*< 

:-"MO"GO 

TO 

4450  \J3=30\GfJ 

TO 

4457 

4450 

IF 

A2$< 

•■"01 "GO 

TO 

4451  \J8"51\G0 

TO 

44  57 

4451 

IF 

A2$< 

"••"LA1  GO 

TO 

4  4  52  \J3-52\G0 

TO 

4457 

4452 

IF 

A2$< 

>"PI"GO 

TO 

4453  \J8~53\G0 

TO 

4457 

4453 

IF 

A2UK. 

>"C1 "80 

TO 

4  4  54  \J8=54\G0 

TO 

4457 

4454 

j ; 

A2$< 

■■"111  "GO 

T  0 

4455  \J8=55\G0 

TO 

4  157 

A  4  5  5 

.1  r- 

A  2  :;>  < 

■•■ill  "GO 

TO 

44  56  \J3~56\G0 

TO 

4457 

4456 

IF 

»'**G0 

ro 

14  57  \J3-57\GG 

TO 

4457 

■157- 


LEDAT.BAS--A  program  to  create  a  binary  vector  file  for 
the  generation  of  characters  and  symbols. 


1005 

OPEN  "R 

.1 

:lei 

FOR 

OUTPUT 

A  S 

FILE  l 

'F1Z(15 

00) 

1010 

READ  B 

1015 

F  =  £+l 

1020 

FOR  1=2 

TO  F 

1030 

REAM  D 

1040 

VF1  (  D-D+f 

1050 

NEXT  I 

1060 

VF1 (1>~I 

< 

1070 

G=VF1 (F 

1075 

G=G41 

1000 

FOR  J  RE- 

■ i 

TO  c 

1090 

READ  E 

1100 

VPKJ)-E 

1110 

NEXT  J 

1120 

READ  H 

1130 

OFKG-i  .1  ) 

sl- 

1140 

FOR  K-Gi 

2 

TO  ( 

3*H+G 

i-2) 

1150 

READ  I 

1160 

vfi  (K  )■--■: 

1170 

NEXT  K 

1180 

CLOSE  ','r 

1 

1200 

END 

1900 

DATA  57 

1905 

DATA  Ifl 

7. 

54  f  8 

3y  104 

•126. 14: 

;:;179y  190 

y  21  7 

1910 

DATA  236 

«2 

52 ,  2 

S2»27i 

:;  „  i)pp  ,  -.? 

.  9  y 

3  4 1 1  3  7 

"iy^OO 

1915 

DATA  440 

« 4 

56  7  4 

■  ;•'!  ,,  ,r1  C 

.  »510»5: 

;3? 

...'  4  2  y  5  b 

•j»571  ti 

02 

1920 

DATA  633 

*6 

46r677r7Qt 

Sy  71  8  k  7: 

'  3  y 

004 r 83 

::«054  »£ 

6  7  ■■' 

077 

1925 

DATA  893 

,9 

1 2  ,  9 

-Qj  ?^  = 

"  y  9  4  3  7  9  <• 

v  --1  v 

°F0y9V 

>>  1000? 

102 

S 

1930 

DATA  104 

7  f 

J.  072 

y 1091 > 

1 1  J.  ?  7 :!. . 

35 

y 1154? 

.  1 7  3 

2000 

DATA  5-1 

•  3 

y  6  ■>   1 

•  3.C-6, 

0  y  -  1  >  2  .» 

ly 

-4»0»0 

>5y-2 

2010 

DATA  12? 

1.  r 

">  F  6  > 

»  4  ..  0  i 

1.2' - 1 ' 

1  ? 

G?-l.i 

,  ...  n  ,  _  1 

2013 

DATA  lr- 

4  , 

0  >  0  r  - 

i  »  0  7  1  y 

2 »  - 1 » 1  * 

07 

-i  y  i .. -: 

2.-1 

2016 

DATA  lr- 

1  , 

3rOn 

>?0 

2020 

DATA  9.0 

,  6 

•  hi 

'-!.-} 

y  .1.  y-4  rC 

f  1 

y  -  1  y  J  y  '. 

-  »  0  7  4  y  1 

7  1  y 

1 

2025 

DATA  1»4 

y  0 

7  1  ,  1 

>  -  1  t  0  y 

0  y  —  5 

2050 

DATA  7..  1 

•  0 

.6?  J 

4  ?  0  y  :l 

*  2  y  --  1  y  3 

>0 

f-A,    1  7 

•2/-1  7  1 

7-4 

y  0  y  0  7  <' 

2060 

DATA  7,0 

,  A 

7  6  ..  1 

-  6  y  0  ? 

1 rOt-6* 

1  y 

6  ;  0  y  0  >  • 

2/3y ly 

.-«j  ? 

0 » 0 1 6  f 

2070 

DATA  6.-0 

•6 

p  6  i  1 

--  6  ;  0  7 

1 ,0,-6, 

0? 

4f3rlt - 

■TvOyO? 

A  v  •- 

~r 

2030 

DATA  11  tOfi 

> ?  5  ,  J 

,  -- 1  ,  :l 

> 1 r~A fC 

7  1 

y -  1  7-1) 

1  ,0.-1 

».!/ 

1  ,  - 1 

2085 

DATA  1?4 

rO 

1  s-  1  , 

1  7  1  7  0 

9 1 » 1 ,-? 

7  0 

y  o>;!.  ••••;• 

2100 

DATA  6*1 

0 

6  •  0  t 

6  7  0  7  1 

y  0  i  -  6  ')  0 

7  0 

-3.  ly-,' 

y  0  y  0  y  6 

•  —  3 

2 1 1  0 

DATA  6fl 

4 

0  r  0  > 

0  y  6  y  1 

y  -  •'!  ?  0  •■  0 

.  2 

f  0  -  1  7  0  , 

-6 «  0  y 2 

,0 

2120 

DATA  6r0 

0 

2  7 1  .« 

2  7  -2  7 

1  f  2  .  0  »  1 

y  »■_ 

.  2  y  1  v  0  » 

4  /  0  y  0  y 

-6 

2130 

DATA  5-1 

0 

6  r  0  t 

6  y  0  y  1 

7  -6y  -3:- 

Oy 

7  >  1  7  1  :•  -1 

v  —  /, 

2 1  A  0 

DATA  3r0» 

0  , 

6  r  1  t 

0  r  -  6  y 

J  y  6  y  0 

2150 

DATA  At  It 

0» 

6  ,  1  r 

3  7  -3  ,- 

I  7  3  :•  3  y  1 

.0 

■-  6 

2160 

DATA  4il. 

0  • 

6  v  1  y 

6y~6y 

1  y  0  y  6  >■  0 

vO 

.  ...x, 

■158- 


2170 

DATA 

2175 

DATA 

2190 

DATA 

2200 

DATA 

2205 

DATA 

2220 

DATA 

i  '.i  "■>  5 

DATA 

2230 

DATA 

2235 

DATA 

2237 

DATA 

2260 

DATA 

2270 

DATA 

2280 

DATA 

2290 

DATA 

2300 

DATA 

2310 

DATA 

2320 

DATA 

2330 

DATA 

2340 

DATA 

23-15 

DATA 

2360 

DATA 

236  5 

DATA 

2380 

DATA 

2390 

DATA 

2395 

DATA 

2410 

DATA 

2415 

DATA 

2430 

DATA 

2440 

DATA 

2445 

DATA 

2447 

DATA 

2470 

DATA 

2475 

DATA 

2490 

DATA 

2495 

DATA 

2510 

DATA 

2520 

DATA 

2530 

DATA 

2540 

DATA 

2550 

DATA 

2560 

DATA 

2570 

DATA 

2580 

DATA 

2590 

DATA 

2600 

DATA 

26.1.0 

DATA 

2620 

DATA 

2630 

DATA 

2635 

DATA 

2640 

DATA 

2650 

DATA 

1  0  r  0  t  1  t  0  i  1  t  - 1  t  I  t  1  t  0  ,  A  t  1  y  1  7  1  f  1  t  4  y  0 
1  y  1  »  -1  y  .1.  y  0  y  -4  y  1  y  -  1  y  "1  y  1  y  -4  y  0  y  0  y  5  y  0 

7tlf0f6tlfAt0rlf2t-lrlr0r - 1 r 1 t -2  r -1 1 1 1 -A t 0 1 0  *  6 » ~3 
llt0f5r0tl*-4t0tlr-ltlrlr0rAflflflrltAf0 

1 1 1 t -1 r 1 y 0  > -4  y 1 y -1 i ~ 1  y  0  y -1 y  2  y 1 v  2  y -2 

8  7  1  «  0  t  6  t  1  ?  4  i  0  y  1  y  2  t  - 1  •  1  -  0  y  -  1  y  1  >  ~2  y  "  .1.  t  1  f  -4  t  0 

0>3»0» 1 »3» -3 

13»0»6f  5»  1  t-1  r  If  1  y-4y0y  1.  y--l  y-1  »  1  t  0  »  - 1 
1  y  1 y  - 1  y  1 »  A  y  0  y  1  y  .1. y  - 1  •  1  »  0  «  - 1  y 1 1  -  1  •  ~  J 
1  »-4f  Orl  »--.!» 1  »0»6»--l 

5  y  0  7  3  ..  0  r  1  y  0  ;  6  «  0  -  -3  y  0  1  1  •  6  y  0  y  0  y  0  »  "6 

7  v  0  v  0  y  6  y  1  y  0  7  -A  t  1  P  2  .«  -2  y  1  f  2  •  0  y  1  y  X  .  2  y  1  y  0  y  A  fOrOt  -6 
A  r  0  7  0  y  6  y  J.  y  3  y  ~6  y   I   y  3  7  6  t  0  t  0  y   "6 

6  y  0  y  0  7  6  :-  J.  y  1  y  "6  y  1  y  2  y  3  r  1  :-  2  y  -3  7  J  y  J   t  6  »  0  f  0  7    V 
4  7  0  y  6  y  6  7  1  7  -6  y  -6  y  0  7  0  »  6  »  1  ?  6  :<  ~6 


6  »  0  ?  0  y  6 


•3  7  1  y  0  .-  "3  7  0  y  0  y  3  y  1  y  3  -  3  >  0  y  0  : 


4  7  0  y  0  7  6  7  1  •••  6  y  0  y  1  y  -6  .»  ~6  ■>  1  7  6  7  0 

5  >  0  >  0  y  4  y 1 y  2  y  2  y 1 y  0  y -6  y  0  y -2 y  0  y  1  •  4  y  0 

I  0  7  0  7  0  7  5  7  1   7  1  7   1  7  I  7  4  7  0  7    1  7  1   7  ~  1   «   1.   7   "'  -  '  "  1   7   1   »  ~  1   t  ~  1 

1  y  -4  y  0  y  1  ,  -  1  y  -  1  y  1  7  0  »  -  2  •  1  •;  6  y  0 


IO»Of  0»6r.l  »6y  0*1 


3yl75.A 


1 


1  f  0 1  -  i  f  1 

4  y  0  y  6  ?  2 1 


■1 


■1  7  .1.  »-4f0»  1  y-lr  1.  fO.'if-l 
Or  1  y6y4y.ly0/~6 


.1  (.?  t  0 1  0  y  1  y  1 .» 1  y  - 1  ?  1  y  4  ?  0  7  .1.  y  1  y  .1.  t  1  1  0  7  2  7  1  y  - 1  »  1 
1  •  ~5  y  0  7  1  7  0  7  2  y  1  y  6  y  0  y  C  y  0  7 -6 
10>0»4»6»lr-4r-4f lfOy-lyl rlr-lt 1 tArO 
1  y  1  y  1  7  1  y  0  7  1  7  1  ?  -  1  t  1  7  1  t  ~A  ,  0  7  C  7  ,:~  ?  -3 

3  y  1  7  6  7  6  7  1  7  -6  7  0  7  0  y  6  7  -6 

1  8  7  0  7  1  /  0  ?  1  7  -  1  y  1  7  1  ■/  0  y  1  y  1  y  1  y  1  j  1  7  4  7  0  y  1  ?  1  ?  1 
irOrl f 1»-1 rlrl»-4»Orl»-3 r-lrli 0»-l ylyly-1 

0  ?  4  7  0  7  1  7  1  y  "  1  y  1  7  0  y  -  J.  -'.  v  "  1  y  - 1  7  1  y  -4  ;  0  y  0  y  5  /  0 

10  y  0  ?  2y  0  y  1  y  4  y  4  y  1  «  0  y  1  y  1  ?  -1  «  1  y  1  •  "4  y  0  y  1  y  -1  »  -1 

1  ?  0  ?  - 1  y  1  y  1  y  - 1  y  1  ■/  4  y  0  »  0  7  J.  y  -  3 

1 1  y  0  y  1  y  0  »  1  y  - 1  y  1  y  1  y  0  y  4  y  1  y  1  ;  1  y  1  v  4  y  0  >  1  y  1  y  ~  1 
1  '  0  7  -4  r  1  y  -1  y  -1  y  1  y  -4  y  0  y  1  y  4  y  6  y  0  ■,  1  y  -6 

5  y  0  y  2  y  Oy 1 y -2  y  2  y 1 y  0  y  2  y 1 y  2  y  2  y  0  y  0 /  -6 

4  *  1  y  2  y  2  y  1  !■  0  y  2  y  1  y  -2  y  2  y  0  .»  2  y  -6 

3  7  0  y  0  y  3  y  1  7  6  7  0  7  0  y  0  »  -3 

5  -  J  v  0  7  1  7  1  7  1  f  0  y  J .  f  0  7  ■-  1  y  1  «  - •  1  y  0  »  0  y  1  r  0 

6  y  0  y  1  y  0  y  1  y  -1  y  0  y  1  ••  0  7  I  ■  1.  -  1  «  C  y  1  ?  0  7  "2  »  0  y  0  7  i. 
5y  0  y  0  y  3  y  .1  y  6y  0  •  0  y  "3  y  -3  r  1  :•  0  -  6  y  0  ?  3  »  "6 
2tlt6t6r0r0i-6 

4  y 1 y  2  f  0  y 1 y  0  y  6  y 1 y -2 • 0  y  0  y  2  y -6 

5  >  0  y  2  y  0  y 1 y -2  y  0  y 1 y  0  y  6  7 1 y  2  y  0  y  0  y  0  y -6 
5  7  0  7  0  y  2  7  1  r  6  y  0  7  0  y  0  7  2  y  1  •  -  6  7  0  y  0  y  6  y  ~  4 
3.  1  76  .3.7  1  y-6y  3y0?6y~6 

3r 0«6y 67 1 7-6.-3; 1 y67-3 

V  7  0  y  3  y  1  7  1  7  -  1  y  1  7  1  y  -  1  y  -  1  :■  1  y  1  1  ~  1  7  1  >   1  7  1  7  j.  7  0  y  2 

.1.  t  J  y  2:  .1.  7  L  yOvOy  1  y~5 

6-0-5  7  2  y 1 y 1 7  0  7 1 y 1 7 ~2 y 1 y 1 y  3  y 1  . 1 t 0  r 0 y 1 y -3 

8  y  0  y  4  y  0  y  1  y  0  y  2  y  1  »  -2  y  0  y  1  r  0  y  -2  y  .1  ;  2  y  0  •  1  y  0  »  "2 


•159- 


2655  DATA  1 , -2y 0» 0 t A t 2 

2660  DA  I  A  6t0rZ'.rA>lrljQrlr2>~At0r 

2670  DATA  ?  t 0 , 1 , 2* 1 r 1 r 1 f 1 t 2 ? 0  y 1 , 1 

2675  DATA  0 t -2 t 0 rltO r 3 t 0 r A t -3 

2680  DATA  5»0»4r0»lv-2»0»lr0»2yl» 

2690  DATA  6  y  Or  2  y  0  r  1  y  0  y  3  y  0  y  0* -1  y  .1  y 

2710  DATA  6f0t2trJrl?0t  -S  >  0  y  0  »  2  »  !.  » 

2720  DATA  ?  y  0  y  2  y  2  y  1  y 2 y 2 y 0 y 0 y -  I  <1  y 

2725  DATA  1  »  2 ,  -2 r 0  r  --1 1 0  y  1  y  0  »  2 »  0  y  3 


Ly2 

1  rO 


l»-2»0*3»0 

-  1  »  1  y  0  »  3 


y  0  »  0  y  2  y  - 

y  0  y  .1.  y  0  y  - 

y  0  v  1  »  0  r  - 

2  j  0  :>  0  y  0  7 


0 1  2  t  0 

0  .  2  y  0 


•160- 


PSTRN.FOR  --A  program  to  send  sequential  data  files  to 
r   the  11  /2p  computer, ........  | 

LOGICAL*  1    FLi ■!  A  f  1  (20?  10)  .■  '."7.7.  •'.  20  ) 

O  JLl  U        „l  U    >'    1  -    1    >    .1    :.) 

TYPE    10 
10  rO!vi"ii'»T(  '    ENTER    FILEi  !AiIE    ■•■    10    A,' 


F.'ir-;n 


r  EF; 

ENTL 


"*  ('  'l 

for;:/;!' 

'.  1  '.•'  i  i 

1  ) 

IFCFLF: 

A  '■'<•'.    i 

r ::  .  i..:' 

50 

cc::ti;: 

j«r 

i  00 

tyte  i 

1  0 

110 

f  ci  :•!■;  at 

(  '  Fr 

IL.  ERA: 

120 

7  V  P  c.  1 

-.-Z-  7  \ 

F"Li!Ai!' 

J.  .•_  >.J 

F- 1)  f\  1 1  A  f 
TYPE  1 

r  1 7  a  .;. ; 

j.  3  0 

;""  !':f,AT 

i  ;  ."; r"  ~  y 

20  0 

?   A  I  '  7 

200 

F  J.'ui m 7 

CIA  4 

II-  '.Mi ■.  ,.:■ 

*  [it!  * 

'  NO  '  ) 

210 

\   ORG  AT 

L.T  I."ir 

ACCEPT 

220 

N  F  i  ■• 

260 


•100 


10    /    {.■'j  ;  ii    J  vO 


r   l :  i  \  i  !  M  I   i.  X  /  ) 

DO    10C0.i-!UH-"l»iN 


L.f-iLL  ir;:;:;; 
Ci-  ...l  a:  :::i 
lT'l;;<u;-<  i:t 


tlG-:Ar.L'F  \  1  ) 

ro   loo-  t-  i ,  err 

Ii"  (  01  G  .  I'-T".  .  A  RUT 
DIu--AI.'UF(  i:< 


'<  0  r  i  -  J.  r 


C« ,";  ,1 


-161- 


1000  COMflilUL 

L-THP     'i:Vt.:' 
20  0  0  ZYOP     'EG  I      LAUADv     F<  I"  C  i ;  Lr  T 

run 


SI.II'.'ROWT  LiJi:     t.'L'l  AC   <[  IAU!    -i>? 
I'll  .ir  i'.'VI  ,',j'     I  ;:(l|    (  ■;  i.v;:  1  ) 

c, '•!.!.    ilall  ( "  i. 6 /'WO/ "  o> 

DG    iOOrJAA.  r  A 

iiXT;-.:o  ■ir-ri:.;;', "  iw.\  ;■ 
i!  (iiXT^o  j;:q,  i  ..:;:;)   c.-g  i  g  :: 
call,   ii-okl*  "  !.>:, 7 ;■"'.: ■:•:?  isur 

GALL    ].  F-'GI  J:.  (  "  x  ..'7  '70  ?  "  1  ) 
Ir  i.  A. -Li  2C -  *  L -C-  .  1  '■      -■    i   ,.•    70 

call   i pes 2:  ( "  1 6 :-"7 ;■•:•?  "o; 


1G< 


1 1 0 


I  AAL 
IPOKIK  "  16; 


EMB 


■162- 


REC.BAS--A  program  to  receieve  sequential  data  from  the  11/34 
computer  and  store  it  in  a  binary  data  file. 


9000  PRINT  "TRANSFERRING  DEC  11/34  DATA  FILES  TO  DEC  11/20  VIRTUAL  DATA 

9010  PRINT  "FILES  (MAXIMUM  VALUE  32767-  MINIMUM  VALUE  0>." 

9015  PRINT  "1023  VALUES  MAXIMUM  TO  BE  TRANSFERRED," 

9020  PRINT  XPK'INT 

9030  PRINT  "HOW  MANY  FILES  TO  RE  TRANSFERRED ■ f \INPUT  NJ 

9050  POR  1=1  TO  Nl 

9060    PRINT    "FILE    I       '  JIJ*     IS'JMNPUT    Lt  <  I ) 

9070    PRINT    L*(I)»"        "  ?  "ON     (Y    OR    N)N\INPUT    Ot 

9080  IF  0$='N'GO  TO  9060  \NEXT  I 

9100  FOR  1  =  3,  TO  Nl 

9105  A*«"RL1 : "SL*<I) S" »DT" 

9110  OPEN  A*  FOR  OUTPUT  AS  FILE  VFJ.%<  1023) 

9120  GOSUB  10000 

9130  CLOSE  VFl 

9140  PRINT  TILE  i       "  J I  ?  '  IS  DONE." 

9200  (4EXT    I 

9250  PRINT  ChR$(7) 

9300  ELL. 

10000  GO  TO  1 00 10 

10010  N-0 

:  .020  DIM  A(15) .B( IS) 

10040  CALL  "  INPB' (B) 

10050  IF  BC1)=1G0  TO  10060 

10055  IF  B(2)-1GQ  TO  10200  \G0  TO  10040 

1 0060  CALL  " DIRA  "  (  1  ? C  , Z ) 

1006--;  VFl  (N)~C 

10065  N=N+1 

10070  A(1)=1\CALL  "PLSB'CA) 

10080  CALL  "INPB' <B) 

100S5  IF  B(1)=1G0  TO  10080 

10090  A<1)*0\CALL  T'LSB"(A) 

10100  GO  TO  10040 

10200  PRINT  N 

10250  N-0 

10300  RETURN 
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OUTPLT.BAS- -A   proaram    to    plot   X-Y    data 


1  DI 

2  Z2 
5  DC-I 
20CO 
2010 
2020 
2030 
2040 

20^;. 

2060 

2110 
212:; 
21^:: 
2140 
2150 
2155 
2156 
2160 
2170 
2210 
2220 

2225 
2230 
n  2  3  p 

2240 
224  2 
2244 
2246 
2  "^  4  7 

O  9  er  g 

2260 
2270 

2275 


2295 

2300 
2305 
2310 
2315 
2320 
2325 
2330 
2340 
2350 


M  B  <  1 5  ) 

=  256 

F  FNA<B9>-B9*1000/4.75 

DIM  X(256) , Y<256) 

P R  .1  NT  '  F I 1  F N  A M !::: '  >  \  I N P U  T  Z 9  * 

OPEN  Z9$  FUR  INPUT  AS  FILE  5  1 

INPUT  tl : ZO* > Z9 i X7 » X6 » Y7  , Y6 

PRINT  Z8*  r  7.9  :-  X7  ,  X6  »  Y 7 »  Y6 

FOR  1=1  TO  Z? 

INPUT  tl :xc I ) , Y< I  ) 

NEXT  I 

CLOSE  41 

A$="I"\GOSUB  4000 

PRINT  *X  SPAN  IS  "  JX7--X6 

PRINT  ' Y  SPAN  IS  " 5 Y7-Y6 

PRINT  "LOWEST  X  VALUE  "  J  X6 

PRINT  ■LOSES:'  Y  MAUL::  "  ?Y6 

PRINT  "NORMALIZE  X  SPAN  TO  HOW 

PRINT  "NORMALIZE  Y  SPAN  TO  HON 

PRINT  "DRAW  AXL'S'fMNPUT  A9$ 

IF  A??:.:-  "Y"GO  TO  2240 

GOSUB  4900  \X1«FNA<1*25*C4)\G0 

X1=0\G03UB  4950  \GOSUB  5000 

GOSUB  4900  \Y1-"FNA(1.25#D4)\G0 


Y1=0\G0SUB  4950  \GOSUB  5000 
PRINT  'OFFSET  X  HOW  MANY  IN CUE 
PRINT  "OFFSET  Y  HOW  MANY  INCHE 
X1-FNA<X1 )\Y1=FNA(Y1 ) 
GOSUB  5000 
C6=C\D6=-D 

"SYMBOL  CHOICE?" 
"AVAILABLE  SYMBOLS  ARE;" 
"SQUARE  <S>  " 
"CIRCLE  (C)" 
° )Rj ANGLE  UP  (Tl  )  " 
"TRIANGLE  DOWN  (T2)  " 
"TRIANGLE  RIGHT  <  T3  >  " 
"TRIANGLE  LEFT''  (  T"4  )  " 
•CROSS  <CR) ' 


MANY  INCHES" 

PMNPUT 

C* 

MANY  INCHES" 

r\ INPUT 

D4 

SUB  5000 

SUB  5000 

3" ?\ INPUT  XI 

S"  y\INF(iT  Yl 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT' 
PRINT'' 
PRINT 
PRINT 
PRINT 
PRINT- 
PRINT 
FOR  T 


"ST  Ah; 
"DOT 
"CHOI 
"WHAT 
1  TO 


(S 
(  D  ) 


)  " 

;\ INPUT  A'ii 
THE  DIMENSION 


\ INPUT  02(1 )\Q2(2)«Q2< 1 ) 


X1=(X<1  )--X.A>*FMACC4>/<X7 
Yl  =  <Y(T)-Y6)*f:NA(D4  )/  (Y7 
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2360  GOSUB  5000 
2370  GOSUB  4000 
2380  NEXT  T 

2385  CALL  '  SI  EP  "  (  -C  •>  30  ,  1  )  \CALL  "  STEP  u  (-D  /  30  »  0  ) 

2386  C=O\D=0 

2390  PRINT  "WOULD  YOU  I..  I K  E  T I C  h  A  R  K  S  •  J  \  I N  P  U  T  0  9  $ 
2400  IF  V9*<>"Y"GG  TO  24  30 

2410  PRINT  "X  VALUE  (-999  FOR  Y  VALUES ) " ? \ INPUT  MS 

2411  IF  M8--999GQ  TO  2430 

2415  Y1=0\X1=  (  MS-X6  )  *FNA  <  C4  )  /  (  X7-X6  )  -!  C6 

2416  GOSUB  5000  \GOSUB  4900  \Xi=C\Yl--20\G0SUB  5000 

2417  GOSUB  4950  \Y1»0\G0SUB  50C0 
24  IS  GO  TO  2410 

2430  PRINT  «Y  VALUE  (-999  FOR  CONTINUE) ' J \INPUT  NO 

2431  IF  N8=-99960  TO  24  70 

2435  X1=0\Y1 = ( NS-Y6 ) *FNA ( B4 ) / ( Y7-Y6 ) +D6 

2436  GOSUB  5000  \GOSUB  4900  \Y1=D\X1=-20\60SUB  5000 

2437  GOSUB  4950  \XJ.-0\GOSUB  5000 
243S  GO  10  2430 

2470  X1=0\Y1=0\G0SUB  5000 

2480  PRIN1  "ANOTHER  PLOT  (Y  OR  N)"J\INPL»T  A9* 

2490  IF  A9*--'Y"G0  TO  2010 

2500  END  


This  section  is  followed  by  the  best-line  routine 
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LASPLT.BAS--A 


program  to  plot  laser  intensity  in  a 
pseudo-three  dimensional  manner. 


900  DIM  B(15) 

1000  PRINT  'X  SIZE"  rMNPUT  X5 

1005  PRINT  "Y  SIZE" J \ INPUT  Y5 

1010  PRINT  "Z  SIZE" J\ INPUT  2! 5 

1015  PRINT  "ANGLE" ?\ INPUT  A 5 

1030  GOSUB  4900 

1035  X1=2£X5\G0SUB  5000 

1040  X1=0\G0SUB  5000 

1045  Y1-2*ZS\60GUE  5000 

1050  Y1=0\G0SUD  5000 

1055  X1~2*X5 

1056  A6~A5*2*3 .141 59/360 

1057  Yl=Xi*SIN<A6)/C0S<A6) 
1060  GOSUB  5000 

1065  GOSUB  4950 

1100  YO-0 

1300  FCK    T8=l  TO  68 

1330  READ  X?Y«Z 

1350  IP  Y=Y060  10  1400 

1360  GOSUB  4950 

1400  Yi--<  Y#SIN<A5*2#3  .  14159/360  )*Y5/1 2)  -fZ* 

14  50  XI  =  <  Y*COS  (  A5*2*3  .  1  4  159/360  )  AY5/12  )  -!  X* 

1460  YO=Y 

1500  GOSUB  5000 

1510  GOSUB  4900 

1520  NEXT  TS 

1530  GOSUB  4950 

1550  X1=0\Y1=0\«OSUE  5000 

1600  END 


Z5/43 
X 15/12 
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